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Abstract 

The majority of the matter in the universe is still unidentihed and under in¬ 
vestigation by both direct and indirect means. Many experiments searching 
for the recoil of dark-matter particles off target nuclei in underground lab¬ 
oratories have established increasingly strong constraints on the mass and 
scattering cross sections of weakly interacting particles, and some have even 
seen hints at a possible signal. Other experiments search for a possible mix¬ 
ing of photons with light scalar or pseudo-scalar particles that could also 
constitute dark matter. Furthermore, annihilation or decay of dark matter 
can contribute to charged cosmic rays, photons at all energies, and neutri¬ 
nos. Many existing and future ground-based and satellite experiments are 
sensitive to such signals. Finally, data from the Large Hadron Collider at 
CERN are scrutinized for missing energy as a signature of new weakly inter¬ 
acting particles that may be related to dark matter. In this review article 
we summarize the status of the held with an emphasis on the complemen¬ 
tarity between direct detection in dedicated laboratory experiments, indirect 
detection in the cosmic radiation, and searches at particle accelerators. 
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6 Towards global analyses 


1. Introduction 


1.1. Objective of this review article 

Only about 5% of the mass-energy content of the universe is composed of 
ordinary baryonic matter, the bulk of which is diffuse gas rather than stars 
and galaxies. Nearly 70% is carried by the so-called dark energy, that is 
mostly observed through its accelerating effect on the expansion of the uni¬ 
verse as observed, e.g., by comparing the luminosity distance and the redshift 
of distant supernovae. The 2011 Nobel Prize in physics was awarded to Saul 
Perlmutter, Brian P. Schmidt, and Adam G. Riess for their groundbreaking 
work that led to the discovery of dark energy. 

In this review we shall discuss the current understanding of the remain¬ 
ing 25% of mass-energy in the universe, the so-called dark matter. Evidence 
for excess gravitational acceleration that cannot be explained by observable 
matter has been found on both small, galactic scales and large, cosmological 
scales. If Newton’s law of gravity, or general relativity, is valid, then the uni¬ 
verse must contain a constituent of unknown nature that betrays its presence 
only though gravitation. This dark matter may or may not be composed of 
particles. If so, the absence of a radiative signal or excess scattering of bary¬ 
onic matter requires that those particles be uncharged and at most weakly 
interacting. Numerous candidate particles have been proposed over the years, 
and a hunt is on to detect them directly via elastic scattering in laboratory 
devices or indirectly through an astronomical decay or annihilation signal. 

Beginning with a historical account, we present the current status of the¬ 
oretical ideas and experimental constraints with an emphasis on complemen¬ 
tarity. Other recent reviews include global accounts of dark matter [DEIE] 
and more focussed publications, e.g., on structure formation |1], baryogenesis 
[3], or indirect searches [HI El El E]- Relevant topics for direct dark-matter 
searches are annual modulation HDl and elastic nuclear recoil [mOSlIIB]. Re¬ 
views of specihc models are available for neutrinos HI, supersymmetry Ha, 
gauginos [161 EH, axinos [18], scalar condensates [m EQ], and asymmetric 
dark matter [21], to only name a few. 

1.2. Observational evidence and alternative explanations 

In 1932, Jan Oort studied the motions of nearby stars and concluded that 
there must be more mass in the local galactic plane than is seen as bright 
stars [22]. Although he interpreted these observations in terms of dim stars, 
they were later taken to be the hrst indication of the presence of dark matter. 
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Figure 1: Rotational velocity curves of seven high-luminosity spiral galaxies of different 
Hubble type. Early-type galaxies have higher peak velocities than later types, but all 
curves remain constant out to large distances (taken from Rubin et al. |26jl. 


Only a year later, Fritz Zwicky made the same claim, based on far better 
data. He assnmed that galaxies in the Coma cluster would be virialized 
and thus found evidence of the total mass of the cluster being more than a 
hundred times that of the stars in the individual galaxies in it [23] • Today’s 
value of the mass excess, based in particular on the now much lower value of 
the Hubble constant of Hq = 67.8 ± 0.9 km/s/Mpc [23] instead of Zwicky’s 
558 in the same units, is about 160 within 5 Mpc |2S|. 

In the 1960s, measurements of the rotation velocity of edge-on spiral 
galaxies became possible. The surprising hnding was that most stars and 
the atomic hydrogen gas have the same orbit velocity irrespective of their 
distance from the center of the galaxy jallEEl E9] and that this extends 
beyond the optical disk. By equating the centrifugal force F = mv‘^/r with 
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the attractive gravitational force F = GnM(< r)m/r‘^, where Gn is Newton’s 
gravitational constant and M(< r) is the enclosed mass, one expects a so- 
called “Keplerian fall-off” of the circular velocity, v oc for tracers 

outside the central luminous disk. The constant circular velocities shown, 
e.g., in Fig. [pimply that the contained mass must increase with radius of the 
galaxy even in the outer regions, where few stars and little gas are found. 
Generally, the total mass is about hve times that seen as baryonic matter. 

Rubin’s observations £t very well to early numerical simulations of galax¬ 
ies composed of 150 to 500 mass points [30]. They demonstrated that cold flat 
disks such as the Milky Way could not survive due to a large-scale, bar-like 
instability beyond a ratio of rotational kinetic energy over potential energy 
of t = 0.14 ± 0.02, unless one postulated a spherical halo with halo-to-disk 
mass ratio between 1 and 2.5. To remain hnite, the halo mass can of course 
at some point no longer increase with r, and the density prohle must fall off 
faster than that of an isothermal sphere, p oc r“^. It is still unknown where 
this happens, and the exact form of the dark-matter density prohle is heavily 
debated. 

The mass distribution in galaxy clusters can also be determined with 
gravitational lensing. So-called strong lensing, which results in multiple im¬ 
ages of a background object, is rare. Weak lensing is based on a statistical 
analysis of the distortions of the shapes of background objects [3I|. The 
consistency of gravitational-lensing results with other types of measurements 
provides convincing evidence that galaxy clusters contain hve times as much 
dark matter than baryonic matter, most of which is in the form of hot, dilute 
gas. 

The most direct evidence for dark matter presented to date comes from 
a system known as the Bullet Cluster [32] • Two individual clusters appear 
to have collided some time ago. Most of the mass, determined using weak 
lensing, is located in the two individual clusters, whereas the intensity distri¬ 
bution of thermal X-ray emission indicates that the baryonic gas is concen¬ 
trated in the collision region, as if it was slowed down when the two clusters 
collided. Hot gas interacts by Coulomb collisions, and so there cannot be an 
unimpeded penetration of the gas haloes of the two clusters. Stars and ap¬ 
parently also dark matter are collisionless and responsive only to collective 
gravitational forces. The dark matter components of the two clusters can 
thus pass through each other. 

Also on cosmological scales one hnds evidence for dark matter. It con¬ 
tributes not only to the total mass-energy content of the universe, but also 
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impacts the growth of density fluctuations. Early in the cosmological evo¬ 
lution, before photons decoupled from baryonic matter, radiation provided 
the bulk of the pressure, but interacted only with ordinary baryonic matter. 
Density fluctuations in the baryonic matter arise from the interplay of grav¬ 
ity and pressure. Dark matter does not react to radiative pressure and thus 
modihes such acoustic oscillations. 

The earliest acoustic fluctuations are measured in the cosmological mi¬ 
crowave background (CMB), whose temperature anisotropies , 6T/T {6, 0) = 
have been measured with unprecedented precision with 
the Wilkinson Microwave Anisotropy Probe (WMAP) [33]. The hrst peak in 
the angular power spectrum of temperature fluctuations is indicative of the 
amount of baryonic matter, whereas the other peaks carry a signature of the 
non-baryonic mass density. More specihcally, the position of the hrst and 
the relative heights of the second and third “acoustic” peaks, indicate the 
universe to be hat, i.e. to have a ratio VL = p/pc of its total density p over the 
“critical value” pc = 3H^/{SnG) very close to unity, and to be dominated by 
cold dark matter (CDM) with a density ratio, Dc, about hve times larger than 
that of ordinary baryonic matter, Db- Data of the subsequent Planck mission, 
when htted with a standard spatially-hat, six-parameter ACDM cosmologi¬ 
cal model, yield CDM and baryonic contents of Qch^ = 0.1186 ± 0.0020 and 
= 0.02226 ± 0.00023, respectively, where h = 0.678 ± 0.009 denotes the 
present-day Hubble constant in units of 100 km/s/Mpc. The baryon density 
inferred from CMB observations is in excellent agreement with the one ob¬ 
tained by comparing measured primordial abundances of light elements with 
the predictions of big-bang nucleosynthesis which at 95% conhdence level is 
given by 0.021 < < 0.025 [Mj. 

After decoupling, baryonic density huctuations can grow, and the uneven 
distribution of galaxies and clusters, known as baryonic acoustic oscillations 
(BAO), is another indicator of dark matter. The Sloan Digital Sky Survey 
(SDSS) performed a survey of 205,443 galaxies at mean redshift of ~ 0.1 
(but also of quasars as distant as z > 5), covering an ehective area of 2417 
square degrees [35]. The three-dimensional data permit studies of BAOs, 
help breaking degeneracies of the two-dimensional CMB, and also provide an 
independent measurement of the total matter content, Qmh = {fib + flc)h = 
0.213 ± 0.023, in agreement with the CMB measurements (note the different 
powers of h). Results from measurements of the CMB and BAO agree in the 
total matter density in the universe being about six times that in baryonic 
matter [36] . 
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Type-la supernova surveys such as the Supernova Legacy Survey (SNLS) 
and that of SDSS provide complementary constraints on the dark-energy 
content of the universe, = ^/(3-ffo)- When interpreted within the stan¬ 
dard flat ACDM cosmology (1 = ffm + their combined data imply 

= 0.295 ± 0.034, which perfectly agrees with the CMB and BAO results 
cited above (note again the different powers of h) [37] . 

It should be noted that the standard ACDM cosmology, involving cold 
dark matter (CDM) and dark energy with similar impact as a cosmological 
constant A, is likewise not without problems. One issue is the missing- 
satellite problem, so named because far fewer dwarf galaxies are observed 
around Milky Way-like galaxies [3S] than are predicted by dark-matter N- 
body simulations. A possible solution could be that only a few dark-matter 
subhaloes experience star formation and evolve to dwarf galaxies on account 
of, e.g., expulsion of baryonic gas by ram pressure [39]. N-body simulations, 
however, point in the opposite direction: most of the massive subhaloes in 
their data are so dense and compact, that they are too big to fail to form stars 
and thus become observable. The velocity dispersion should be signihcantly 
larger than that observed in real dwarf galaxies. 

Variations of the laws of gravity have been proposed as a means to explain 
all observational hndings without dark matter. Developed in 1983, modihed 
Newtonian dynamics (MOND) posits that below a certain level of gravita¬ 
tional acceleration, oq — 2 ■ 10“^° ms“^, the gravitational force is no longer 
proportional to the acceleration, but the acceleration squared BDl. Origi¬ 
nally a purely phenomenological nonrelativistic model designed to reproduce 
flat rotation curves, a relativistic formulation of MOND, called TeVeS, is now 
available [H] , that also predicts gravitational lensing and furthermore evades 
the problem of superluminal scalar-wave propagation of previous attempts. 
TeVeS is built on the notion that the dynamics of matter helds is governed 
by a metric tensor — 2UfjUy sinh(20) that is different from the 

metric appearing in the Einstein-Hilbert action. In particular, a timelike 
vector held satisfying g^^'^u^Uy = — 1 and a scalar held 0 were introduced 
to exponentially stretch and shrink the Einstein metric in the spacetime di¬ 
rections orthogonal and parallel to respectively. The scalar held action 
is governed by a dimensionless function and a free length scale, that can be 
related to the parameters /i and oq of MOND in the quasistatic limit. In 
addition, the theory contains a scalar and a vector coupling constant as free 
parameters. 

Although most physicists prefer dark matter over MOND as an explana- 


tion of the excess gravitational acceleration in galaxies and the universe, the 
debate continues 021 . As an example, objections based on the rapid insta¬ 
bility of spherical systems such as stars |13] were alleviated by a generalized 
version of TeVeS |15 . On the other hand, it is worth noting that the Bullet 
Cluster provides evidence for dark matter independent of the particulars of 
the laws of gravity |1S]. MOND, or TeVeS, therefore, certainly has trouble 
accounting for the Bullet Cluster, but its predictions must nevertheless be 
tested, for example by measuring the orbits of binary stars of large separa¬ 
tion, for which the gravitational acceleration is below oq. 

1.3. Fundamental requirements for dark matter particles 

Let us suppose that dark matter is composed of individual particles with 
rest mass mo- For hot dark matter, mo (? is less than the energy scale of the 
universe shortly before recombination, and the particle was relativistic while 
structures started to grow, suppressing growth on small scales by smearing 
out such density fluctuations. Warm dark matter has a mass, and hence a 
velocity, just large enough that structure formation on the scale of a proto¬ 
galaxy is possible. Only cold dark matter is massive enough to permit growth 
of small-scale density fluctuations. 

Simulations can be used to follow the evolution of structures after recom¬ 
bination, when radiative pressure is negligible. Structure formation in this 
phase is hierarchical, meaning that small structures form hrst and larger ones 
later. The matter distribution at redshift z = 0 still carries imprints of the 
fluctuation spectrum at the time of recombination. Cold dark matter seems 
to £t best, and hot dark matter as dominant factor is virtually excluded |16]. 
Measurements of the Lyman-a forest, small-scale gas structures at redshifts 
of a few, are incompatible with a particle mass below the keV scale HU- 

Very massive cold objects, for example planets, brown dwarfs, or primor¬ 
dial black holes, which are colloquially subsumed under the name Massive 
Compact Halo Objects (MACHOs), have been extensively searched for with 
microlensing experiments and found too few in number to account for dark 
matter [IHl SSI ED] • 

Dark matter cannot be baryonic, because that would be in conflict with 
the yield of big-bang nucleosynthesis IS] and results of a CMB fluctuation 
analysis. The mass range for which primordial black holes can constitute dark 
matter is also quite restricted to roughly the mass range between 10^’^ Mq 
and 10^^ Mq [52] • This makes it very likely that dark matter is a new type 
of particle. 
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A consequence of non-baryonicity is color neutrality, meaning dark matter 
particles may not engage in strong interactions, otherwise that particle would 
just act as baryonic matter. In its lightest form, it has to be electrically 
neutral, lest it be easily visible instead of being dark. Note that this does 
not preclude the existence of heavier charged siblings that may have existed 
in the early universe. 

Weak interactions are not excluded. A coupling to the electroweak gauge 
bosons and Z must be somewhat weaker than that of other Standard- 
Model particles, though, otherwise direct-detection experiments would al¬ 
ready have seen such interactions. 

The stability of at least the lightest dark-matter particles must be high or, 
conversely, the self-interactions rather weak. The lack of impact on the dark- 
matter distribution of cosmic events that do influence the baryonic-matter 
distribution, such as the merger of the Bullet Cluster imply that the dark- 
matter particles have to be highly collisionless, much more so than the intra¬ 
cluster gas that is itself collisionless in the usual dehnition of a preponderance 
of collective interactions over two-body collisions. In fact, these observations 
have been used to put upper limits on the dark matter self-interaction scat¬ 
tering cross section per dark-matter mass mx due to the resulting drag force 
per mass unit, F^/mx ~ {o's/mx)v‘^PX) with v and px the relative velocity 
and mass density, respectively. These constraints are currently of the order 
o's/^x ^ Icm^/g [S3]. Interestingly, a signihcant offset between the cen¬ 
ter of the dark-matter halo and the star distribution of one of the central 
galaxies in the galaxy cluster Abell 3827 has been claimed and interpreted in 
terms of a dark-matter self-interaction scattering cross section of the order 
of as/mx 10 "^cm^/g [SI], but this is currently controversial |55j . 

Depending on the nature of the decay products, astrophysical lower limits 
on the decay timescale are > 10^® s [e.g. [56]. Self-annihilations must also 
be rare. Assuming that the particles were in thermodynamic equilibrium 
very early in the universe, the present-day average density of dark matter 
requires some loss in particle number on account of, e.g., annihilation. If S- 
wave processes dominate, i.e. the product of velocity and annihilation cross 
section, (cxxv), is constant, then the X—particle mass density in units of 
the critical density today is given by 


1/2 

9f 
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where 10 < < 100 is the number of relativistic degrees of freedom at the 

freeze-out temperature, Tj ~ mx/20, below which the dark-matter particles 
of mass mx cease to be in thermal equilibrium, as can be shown by solving 
the appropriate Boltzmann equation. Furthermore, Hq and Tq are the Hubble 
rate and CMB temperature today and Mpi = 1.22 x 10^® GeV is the Planck 
mass. This thermal freeze-out process and its scaling with the annihilation 
cross section can qualitatively be understood as follows: The rate of change of 
the dark-matter number density, nx, due to expansion of the universe equals 
—3H{T)nx with H{T) the Hubble rate at temperature T, whereas the one 
due to annihilation and creation from inverse processes is proportional to 
{o'xxv)n\. For T > mx the number density, nx, is given by the thermal 
equilibrium abundance, and for T -C mx annihilation and inverse processes 
become inefficient so that nx oc due to the expansion of the universe. 
In the transition region T ~ mx both rates become comparable. Equating 
these rates and using the Friedmann equation for H{T) then leads to the 
scaling in Eq. Q. This also implies that the annihilation cross section has 
to be equal to or larger than the thermal cross section [571 EH] if the particle 
is to constitute all or part of the dark matter, respectively. 

In principle, the particles constituting CDM could have any spin. If they 
happen to be fermionic and non-degenerate, an interesting lower bound on 
the fermion mass, known as the Tremaine-Gunn bound [3^1, can be derived. 
Fermions satisfy Pauli’s principle, and so the fermionic occupation numbers 
are always below the Boltzmann limit, feq{E) < g exp[{fi — E)/T]. With 
fi = m one hnds 


p ~ mn < gm 



~ gm 


f m^al 
\ bvr 


3/2 


( 2 ) 


where in the last step we have assumed a Gaussian distribution to compute 
the velocity dispersion cr^ in one direction, ^m{v^) ~ ^ma^ = T/2. We 
know that individual dark matter haloes exist. In the simplest form they 
are isothermal, = const., for which the hydrostatic equilibrium implies 
p oc r“^. Inserting these expressions one hnds 


m 


> 


(27r)i/8 


{gG^a„r‘^yi‘^ 


~ 1.5 (yf 


- 1/4 


/1000 km s /Mpc\ 
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(T7, 




eV. 


(3) 


Some variation in the pre-factors is to be expected if not all assumptions 
are met, but Eq. (|^ is accurate to within a factor of a few. This limit is 
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most severe for small haloes with low velocity dispersion. For dwarf satellite 
galaxies of the Milky Way we may estimate r 500 pc and ~ 50 km/s, 
resulting in m > (150 eV). Obviously, very light fermionic particles 

constituting hot or warm dark matter appear inconsistent with the existence 
of small dark matter haloes. 

The Standard Model of particle physics does not contain any particle 
that meets all these requirements. In the following we list three popular 
types of extensions to the Standard Model in which dark-matter candidates 
are discussed. An overview over these and other candidates is given in fig. 

1. In supersymmetric extensions of the Standard Model the lightest super- 
symmetric particle (LSP) is stable on account of a specihc symmetry. 
Such particles typically have masses above a few GeV and are thus 
candidates for CDM. Among these so-called weakly interacting mas¬ 
sive particles (WIMPs), one hnds the supersymmetric partners of the 
gauge bosons and the Higgs boson, i.e. a neutralino, of a neutrino, i.e. 
a sneutrino, or even of the graviton, i.e. the gravitino. 

2. Also subsumed as WIMPS are Kaluza-Klein excitations found in higher¬ 
dimensional extensions. Again, the lightest particle of this type can be 
a dark-matter candidate. 

3. Axion-like Particles (ALPs) are generally very light, with masses below 
the MeV scale and typically much less, down to ~ 10“® eV. In contrast 
to WIMPs, they are not produced through thermal freeze-out, but form 
a Bose-Einstein condensate with very high occupation numbers. 

Of particular interest and subject to controversies is the distribution func¬ 
tion of dark-matter particles. This is not an academic question, because the 
velocity distribution of particle in the barycenter of the solar system is the 
decisive input parameter for the interpretation of direct-detection experi¬ 
ments. The density prohles of dark-matter haloes and clumping on small 
scales determine the yield of dark-matter annihilation for indirect-detection 
experiments [e.g. |6T], as the rate depends on the square of the density, p. 
For any fluctuating quantity the volume average of the density squared is 
amplified. 


P = Po + Sp ^ (P^) = (pD + {{Spr) > (pD . (4) 


Simulations suggest that the mass distribution of dark matter clumps is well 
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Figure 2: Summary of dark matter 
candidates in double-logarithmic dis¬ 
play of mass mx and typical scatter¬ 
ing cross section with ordinary matter. 
Hot, warm, and cold dark matter are 
indicated in red, pink and blue, respec¬ 
tively. The gravitino is the supersym¬ 
metric partner of the graviton, whereas 
neutralinos are candidates for the LSP. 
ADM denotes asymmetric dark matter, 
to be briefly discussed in Sect. |2.4[ and 
ordinary as well as sterile neutrinos are 
candidates for hot and warm dark mat¬ 
ter. Finally, SIMP stands for strongly 
interacting massive particles and WIM- 
PZILLAS consitute supermassive dark 
matter that can only be created at 
GUT energy scales. Taken from Baer 
et al. [50] . 


represented by a power law 


dn{M) 

dM 


TIqM 


(5) 


whose index is approximately / = 2 between 10“® Mq and 10 Mq [02], and 
somewhat flatter, / ~ 1.8, around 10® Mq [e.g. EB]- We note that the smallest 
clump mass is determined by kinetic decoupling of the WIMPs and thus by 
their microphysics and can range between 10 Mq and ~ 10 ^ Mq [M] . 
Apart from the low-mass cut-off of the dark-matter clump distribution, the 
total number density of clumps also depends on the shape of the mass func¬ 
tion which in turn depends on the poorly known survival probability of such 
microhaloes [e.g. Eg. 

For a dark matter clump of characteristic size, r^, and density, the 
annihilation rate can be written as [Ml EZ] 


Qaiin CC CC M' 


( 6 ) 


For a cuspy density profile ps should also be proportional to the mass of the 
clump, M. Then, Qann may increase faster than linearly with mass [Ml ES], 
i.e. d > 1, and the annihilation rate per mass interval peaks at the high-mass 
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end of the distribution, 


dQ 

dlnM 


QcM 


dn{M) 

dM 


oc . 


(7) 


If the power-law index in equation ([^ is negative, d -1- 1 — / < 0, then the 
dark matter annihilation is dominated by the numerous low-mass clouds, and 
we can treat clumping with a simple boosting factor. On the other hand, if 
d -|- 1 — / > 0, then the few massive clumps dominate electron production, 
and the spatial distribution of annihilation events will no longer reflect the 
global dark-matter density prohle. 

On large scales, the density prohles of dark-matter haloes are likewise not 
well known. Whereas various mathematical expressions can be adapted to 
£t the density prohles in the outer regions of the haloes, the behavior in the 
central regions is subject to a controversy known as the cusp/core problem 
m- The rotation curves of the innermost regions of dwarf galaxies suggest 
that the density prohle levels oh and forms a core. Observations of galaxy 
clusters likewise indicate a shallow mass prohle. If the rotation velocity falls 
oh faster than r for r —)■ 0, then the mass density, p(r), cannot increase for 
r —)■ 0. Such a prohle is well approximated as 


Pisoir) 


Ps 



( 8 ) 


which turns into an isothermal prohle for r > Vg. In contrast, N-body dark 
matter simulations predict dark matter prohles that behave as p{r) oc r~^ 
or as p{r) oc r~^'^ for small r. The most widely used density prohle, the 
Navarro-Frenk-While (NFW) prohle [71], is an adaption to simulation results 
and given by 

Pnfw(’") = — . ^ '2 ’ (^) 


1 + ^ 
rs 


where r* is a scale radius and pg a characteristic density. More recent simula¬ 
tions suggest a slight modihcation in the innermost regions, which is known 
as an Einasto density prohle [72] . 


PEinasto (j"') 


Ps exp 



exp 




with a ~ 0.16 . 


( 10 ) 
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The slope of this prohle depends on the radius r as 


din p{r) 


d\nr 



( 11 ) 


Warm dark matter would lead to a flattening of the central density prohles of 
dark-matter haloes compared to the ones discussed above for CDM. A class of 
explanations for the cusp/core problem focusses on the neglect of baryons in 
dark-matter simulations [73]. Star formation and supernova explosions may 
remove low-angular-momentum gas m- Alternatively, angular-momentum 
transfer between baryons and dark matter may lead to heating of dark mat¬ 
ter. One also needs to verify that residual collisionality in N-body simulations 
does not artihcially shape the central density prohle [7511761I77]. 

2. Candidates 

2.1. Neutrinos 

Standard-Model neutrinos can not constitute a signihcant part of dark 
matter because current upper limits on their mass imply that they would 
freely stream on scales of many Mpc and hence wash out the density huc- 
tuations observed at these scales. Many extensions of the Standard Model 
contain one neutrino state per lepton generation which are singlets under the 
electroweak gauge group. Such states are generally termed sterile neutrinos^ 
although they are not sterile in the strict sense because they may be sub¬ 
ject to other, “hidden” interactions, and the mass eigenstates are in general 
mixtures of these electroweak singlet states and the active neutrinos. If the 
mass of the heaviest neutrino state is on the order of the Grand-Unihcation 
scale, then one can show that the diagonalization of the 6x6 mass matrix 
of active and sterile neutrinos can naturally provide the small, sub-eV mass 
scale of the known active neutrinos. This concept is known as the seesaw 
mechanism. 

If the lightest sterile neutrino has a keV mass scale, it is a warm dark- 
matter candidate na. The couplings of sterile neutrinos are usually too 
small to have been in thermal equilibrium in the early universe. However, 
if they mix with active neutrinos, they can be produced through oscillations 
with active neutrinos which in turn have been in thermal equilibrium at tem¬ 
peratures above an MeV in the early universe. If the total lepton number in 
the neutrino sector vanishes, this process is known as the Dodelson-Widrow 
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Figure 3: Summary of current constraints on sterile-neutrino mass and mixing angle with 
active neutrinos, from Ref. |81j . Red squares mark the values implied if the putative 
3.55-keV line were due to the decay of a sterile neutrino, constituting the entire dark 
matter, into an active neutrino and a photon through a one-loop diagram involving a 
W~. Left panel: The colored area on the upper right is ruled out by various observations 
of diffuse X—ray emissions. The red hatched area marked “Pulsar Kicks” would give the 
observed kicks to newly-born neutron stars due to asymmetric emission of sterile neutrinos. 
Red lines denote parameter combinations, for various lepton numbers as indicated, for 
which sterile neutrinos would be produced from active neutrinos through oscillations and 
would constitute all dark matter, known as the Dodelson-Widrow scenario [T^l. The 
general lower bound on fermion dark-matter mass marked “Tremaine-Gunn” is discussed 
in Sect.[^ see Eq. @ and also Ref. [82] for more information. Right panel: Recent Jf—ray 
constraints based on observations of M31 by Chandra and XMM-Newton and of Ursa 
Minor by Suzaku. The Dodelson-Widrow band is the same as on the left panel (L = 0). 
See also Ref. [13| for more information. 


scenario [TU] , whereas the case of non-vanishing lepton number is often called 
the Shi-Fuller scenario [20] • One-loop diagrams involving a W~ give rise to 
decays into an active neutrino and a photon. If sterile neutrinos signihcantly 
contribute to dark matter, these decays can lead to detectable diffuse X—ray 
emission, measurements of which rule out a considerable part of the param¬ 
eter space of sterile neutrinos. Figure [^summarizes constraints on the mass, 
nis, and the sterile-active mixing angle, 0. Also shown in Fig. [^are the im¬ 
plications of interpreting in the context of sterile-neutrino dark matter the 
recent putative observation of a 3.55-keV line. 
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2.2. Weakly interacting sub-eV (or slim) particles (WISPs) 

Many extensions of the Standard Model of particle physics, in particnlar 
scenarios based on snpergravity or snperstrings, predict a hidden sector of 
new particles interacting only very weakly with Standard-Model particles. 
Snch scenarios do not necessarily only contain Weakly Interacting Massive 
Particles (WIMPs), new heavy states at the TeV scale and above, some 
of which are candidates for the dark matter, bnt often also predict Weakly 
Interacting Snb-eV (or Slim) Particles (WISPs) that can conple to the photon 
field, [HI]. Well-known examples are psendo-scalar axions and axion-like 
particles (ALPs), a, and hidden photons that kinetically mix with photons. 
At the high end of the mass spectrum, various particle-physics models also 
predict the existence of non-elementary particle states that can be either one¬ 
dimensional topological defects such as monopoles or non-topological solitons 
such as condensations of bosonic states and so-called Q-balls. In the following 
we briefly discuss the motivation for and constraints on some of the lighter 
dark-matter candidates. 


2.2.1. Axion-like particles (ALPs) 

In addition to the standard QCD Lagrange density, QCD allows a term 
violating CP symmetry that is of the form 


^0 - 


Act 5 


( 12 ) 


where 6*“^ is the gluonic field-strength tensor with a = 1, • • •, 8, and as 
its dual, and 6* is a dimensionless number whose size has to be experimentally 
determined. Experimental limits on the CP symmetry-violating electric- 
dipole moment of the neutron requires the parameter 9 = 9 + arg det M to 
be smaller than ~ 10“^°, where M is the quark mass matrix. There is a 
priori no theoretical reason why this number should not be of order unity, 
and so this small value is known as the strong CP problem. 

A possible solution to this problem is to promote 0 to a dynamical pseudo- 
scalar axion field., a, whose expectation value can be dynamically driven to 
zero. This behavior can be achieved, if the low-energy effective Lagrange 
density for the axion has the form 


Ca = -d^ad^a + 


+ 8 ^ 




2 „2 


(13) 
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where s is a model-dependent number, fa is the Peccei-Quinn energy scale 
and is the electromagnetic-held tensor with its dual Fxa/2. 

If the effective axion mass, nia, vanishes, the effective action based on the La¬ 
grange density (cf. Eq. 13) is invariant under the Peccei-Quinn shift symme¬ 
try, a —)■ a-|-const. This invariance follows from the fact that the parity-odd 
terms involving held strengths in Eq. (13) can be written as divergence of a 
chiral fermion current. The divergence can then be transferred to the axion 
held by partial integration, and only derivatives of the axion held remain. 
Furthermore, by comparing Eqs. (12) and (13) one sees that 9' = 9 -\- sa/fa 
plays the role of an angle, that in the early universe within one causal volume 
will take a random value of order unity. This aspect will play a role for the 


axion relic density that we will discuss in Sect. 2.2.3 
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Figure 4: Constrained and still testable regions of photon-ALP mixing in the ma — ga-y 
plane. Experimentally excluded regions are shown in black, constraints from astronomical 
observations in gray, and constraints from astrophysical and cosmological arguments in 
blue; the sensitivity of planned experiments is indicated in light green. Plotted in red are 
the parameter regions in which ALPs could account for all dark matter, see Sect. |2.2.3| 
The KSVZ axion is one of the original QCD axion models which roughly exhibit the scaling 
in Eq. From Ref. |86) . 
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The coupling of Axion-Like Particles (ALPs) to photons can now be put 
into the form 

rv... ~ rv... 

aE-B, (14) 


p _ ^em 






2vr/a 


where in the last expression we have rephrased F^yF^'^ in terms of the elec¬ 
tric and magnetic held strengths E and B, respectively, and the often-used 
convention 

^ (15) 


9a'y — 


2vr/„ 


1 


was applied. Eq. (14) is the basis to derive many oscillation phenomena 


between photons and ALPs, as we shall see in Sect. 4.6 


The ALPs couplings to gluons give rise to mixing between ALPs and 
neutral pions, tt^. In the original QCD axion scenario, the chiral symmetry 
breaking of the strong interactions which gives rise to the neutral-pion mass, 
also leads to an axion potential which at zero temperature has the form 


m 


r(o) AtcD ’ 


( 16 ) 


where Aqcd — 215 MeV is the conhnement scale. Expansion for a//a <C 1 
then yields an axion mass 


m. 


~ 0.6 


10iOGeV\ 

- 7 - I meV , 

Ja J 


(17) 


which is thus inversely proportional to the Peccei-Quinn scale. We mention in 
passing that in so-called natural inflation scenarios approximately massless 
(pseudo) Nambu-Goldstone bosons such as the axion can also play the role 


of the inflaton with a potential of the form Eq. (16), with Aqcd replaced by 
a more general energy scale. 

Fig-i shows a compilation of current constraints on and sensitivities of 
planned experiments to photon-ALP mixing in the rua — Qa'y plane. 

In the laboratory photon-ALPs mixing can be probed with so-called shin¬ 
ing light through a wall experiments. Photons of a laser beam are partly 
converted to ALPs in a strong magnetic held in front of a wall and then 
reconverted behind the wall by a similar magnetic held in an optical cavity 
of high Q value. One such experiment, the Axion-Like Particle Search^ is 


^Alternatively called Any Light Particle Search (ALPS), https://alps.desy.de 
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operated at DESY and uses a 5-T magnetic field and an optical cavity of 
8.4m length. 


2.2.2. Hidden photons and other WISPs 

A hidden-photon field, describes a hidden f/(l)-symmetry group and 
mixes with the photon through a Lagrange density of the form 

1 1 sin V cos^ V 

Cx, = 

( 18 ) 

where is the hidden-photon field-strength tensor, m.^/ the hidden-photon 
mass, and y a dimensionless mixing parameter; is the electromagnetic 
current. Typical values for the mixing parameter range from ~ 10“^ down 
to 10-1^ 
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Figure 5: Constrained and still testable regions of photon/hidden-photon mixing in the 
my — X plane. Experimentally excluded regions are shown in green and gray, the other 
colors indicate constraints from astronomical and astrophysical observations, as indicated. 
Constraints from CMB distortions are shown in black. Below the thin red line the hidden 
photon would account for all the dark matter observed by non-thermal production in the 
early universe. From Ref. [55] . 

Fig-i shows a compilation of current constraints on, and sensitivities of, 
planned experiments to photon/hidden-photon mixing in the my — x plane. 
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The ALPS experiment can also be used to probe photon/hidden-photon 
mixing when operated without magnetic held. The resulting constraints are 
included in Fig. 


2.2.3. WISPs as dark matter 

WISPs are typically very weakly coupled and very light. If they have 
ever been in thermal equilibrium in the early universe, they would have 
thermally decoupled already at very high temperatures. In this case, as long 
as they have not formed a Bose-Einstein condensate, they would contribute 
to hot dark matter, and their mass density today would be far too small 
to explain dark matter. Therefore, if WISPs signihcantly contribute to dark 
matter today, very likely they have been produced non-thermally in the early 
universe and/or formed a Bose-Einstein condensate. 

Let us discuss non-thermal WISP production modes on the specihc ex¬ 
ample of the QCD axion. There are essentially two potential contributions 
to non-thermal relic axions, and both are related to breaking of the Peccei- 
Quinn symmetry: First, once this symmetry is broken, there will be a po¬ 
tential energy whose zero-temperature limit we encountered in Eq. (16). In 


general, the axion held, a(f,r), will not be at the minimum of this poten¬ 
tial, resulting in oscillations around this minimum associated with kinetic 
and potential energy. Oscillations of the homogeneous mode of scalar or 
pseudo-scalar helds manifest themselves as non-relativistic dark matter with 
an initial energy density given by 


PaAT) ^ ^ Olo < 


(19) 


where ma{T) is the temperature-dependent axion mass, oq is the value of the 
axion held at the time of symmetry breaking, s is a dimensionless parameter 
that depends on the axion model, and Oafl = sa^/fa the corresponding angle. 
Since the resulting energy density is proportional to the square of the devi¬ 
ation from the equilibrium position, 6a,o, fbis process is also known as the 
axion misalignment mechanism. The resulting energy density today can be 
obtained by solving an equation of motion analogous to that for the inhaton 
held, 

d + 2,H{t)a + m?Jfr{t)]a = Q , (20) 

which is a damped-oscillator equation with a time-dependent mass ma[T(f)]. 


21 




In terms of the critical energy density, one obtains 


~ 0.11 


a 


fa 


5 X IQii OI q GeV 


1.184 1.184 

,2 ( 12 /leV 


~ 0.110 


a,0 


rrir, 


( 21 ) 


where in the second expression we have used the relation Eq. 0 between 
the zero-temperature mass, nia, and fa for QCD axions. A generalization of 
Eq. (1^ to ALPs gives 




1/2 


fa 


1011 GeV 


2 /n \ 2 
^a,0 


TT 


( 22 ) 


If the ALP or axion held becomes massive before inhation, the initial mis¬ 
alignment angle, 6a,o, is the same as that in the visible universe today, whereas 
if the held becomes massive only after the end of inhation, the misalignment 
angle will be diherent in diherent Hubble volumes at the end of inhation and 
must be averaged to obtain today’s energy density, {OI q) = vr^/S. 

A second contribution to relic axions arises because PQ-symmetry break¬ 
ing also gives rise to axionic strings, since breaking of a U{1) symmetry al¬ 
ways gives rise to strings. The Higgs-Kibble mechanism predicts that about 
one string will be produced per causal volume. If PQ-symmetry breaking 
occurs after inhation, such that the strings are not inhated away, the decay 
of these strings will then give rise to the contribution 


nih^ ~ 0.11 


/ 400 fieV \ 
\ ma J 


(23) 


It is interesting to note that if the inhation scale is close to the upper 
limit given by the tensor-to-scalar ratio, then most likely the ALP scale 
/a ^ hfi — 10^^ GeV because otherwise axion isocurvature huctuations would 
be produced during inhation, inconsistent with upper limits from GMB ob¬ 
servations. From Eq. (21) with (0a,o/7r)^ —t 1/3 and Eq. (23) one would 
deduce a rather narrow window of 10® GeV < A < 10^^ GeV in which ALPs 
or axions could represent the bulk of dark matter. According to Eq. (17) the 
axion mass would then be in the range 10“® eV < rua ^ 10“® eV, whereas 
ALP masses could be in the eV range. 

If WISPs make up a signihcant part of dark matter, they can also be 
experimentally detected with so-called haloscopes which typically consist of 
a resonant cavity that can be tuned to the photon energy corresponding 
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to the WISP mass. The leading Axion Dark Matter eXperimenlj^ (ADMX) 
uses a cryogenically cooled high-Q value tunable microwave cavity that is 
immersed in a magnetic held of about 8 T. The resulting constraints on 
ALP parameters are shown in Fig. An extension to higher frequencies, 
ADMX-HF, is planned. 

Very recent ideas focus on dish antennas to detect the conversion into 
radio photons of WISP dark matter in the mass interval 10“® eV < mwisp ^ 
leV [87]. 

Finally, due to their, albeit very weak, coupling to photons and, possibly, 
other Standard-Model particles, WISPs can also provide potential indirect 
signatures via their decays in Standard-Model particles. For example, the 


ALP coupling to two photons (cf. Eq. 14) leads to a lifetime 


Ta = 


bdvr 


~ 5.1 X 10 


26 


fa 


1010 GeV 


evy 

maj 


s. 


(24) 


For the preferred parameters discussed above this is obviously much larger 
than the age of the universe. If ALPs constitute a fraction of dark matter, for 
rua ~ 7.1 keV and lOH GeV < /a < 10^® GeV, the lifetime may be sufficiently 
short to explain, for example, the recent indications for a 3.55-keV photon 
line that we discuss in Sect. 14.21 


2.3. Weakly interacting massive particles (WIMPs) 

Supersymmetry (SUSY) relates fermions to bosons whose spins differ by 
1/2, and they appear as so-called super-partners in super-multiplets. Fur¬ 
thermore, SUSY commutes with the energy-momentum operator. Among 
other things this implies that, for unbroken supersymmetry, the rest masses 
of all members of a super-multiplet must be equal. Since such pairings are 
obviously not realized in nature, SUSY must be broken. An important fea¬ 
ture, and one of its main motivations, is that even broken SUSY stabilizes 
scalar masses which would otherwise be unstable to radiative corrections, be¬ 
cause SUSY links the scalar mass to the mass of the corresponding fermionic 
partner which is protected by gauge symmetry. Furthermore, supersymmet¬ 
ric extensions of the Standard Model typically have an additional discrete 
R-symmetry, also called R-parity. Since supersymmetric partners of ordi¬ 
nary Standard Model particles are odd, the lightest SUSY particle, known as 


^http://www.phys.Washington.edu/groups/admx/home.html 
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the lightest supersymmetric particle (LSP), is thus stable against decay and 
constitutes a candidate for dark matter, called a weakly interacting massive 
particle (WIMP). In particular, the super-partners of the neutral electroweak 
gauge bosons, the photon and and of the two neutral Higgs-boson states 
present in SUSY models, or linear combinations of these four states, are 
promising WIMP candidates because of their weak interactions. They are 
known as neutralinos. 

WIMP candidates also occur in extensions of the Standard Model contain¬ 
ing compactified extra dimensions. Such scenarios are motivated by string- 
theory models which require extra dimensions to be anomaly-free, and also 
by the fact that in a flat space of n extra dimensions and a volume 14 the 
fundamental, higher-dimensional gravity scale, M*, is related to the Planck 
mass, Mpi, observed in fonr dimensions by the relation 

= (25) 

This allows to “tnne” 14 snch that the fnndamental gravity scale, M*, can be 
aronnd a TeV, thus avoiding the hnge and nnexplained difference between the 
electroweak scale and gravity. Assnming that the wave fnnction factorizes 
into an ordinary three-dimensional part and a part that only depends on the 
coordinates of the extra dimensions, Sj, the latter must be single-valued and 
thus must have the form 


exp 



2'KniSi \ 
Tn ) 


1 


with Hi G Z. The dispersion relation of a particle with momentnm p 3 in the 
ordinary (inhnitely extended) three spatial dimensions then reads 


B"(p) 



(26) 


From the point of view of ordinary three-dimensional space, there may be so- 
called Kaluza-Klein excitations of the qnantnm helds describing elementary 
particles, with mass contribntions equal to multiples of 27r/r„, 


27r n 

— ~ 6 X 10^ 
Vr, 




GeV. 


(27) 
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If the compact space of extra dimensions has a reflection symmetry, as is the 
case in the example above, the so-called Kaluza-Klein parity, (—is 
conserved at tree level. The lightest Kaluza-Klein state is thus stable or long 
lived which makes it a dark-matter candidate. 

2.4- Alternatives 

The abundance of WIMPs discussed above is governed by the thermal 
freeze-out paradigm, in which WIMPs drop out of thermal equilibrium at 
temperatures T < mx/20, as we have discussed in Sect. However, there 
are also dark-matter candidates which have never been in thermal equilib¬ 
rium in the early universe. This is typically the case for sterile neutrinos 
and WISPs. But whereas sterile neutrinos and WISPs tend to be relatively 
light and are produced by mixing with Standard-Model particles or through 
the dynamics of a scalar or pseudo-scalar held, so-called feebly interacting 
massive particles (FIMPs) with mass comparable to the electroweak scale, 
alternatively also called frozen-in massive particles, could be produced by 
collisions or decays of Standard-Model particles. If such particles start from 
a vanishing initial abundance, they are said to be produced by freeze-in, 
in contrast to the freeze-out scenario discussed above, in which dark-matter 
particles have interacted frequently enough with Standard-Model particles to 
have initially been in thermal equilibrium. Let us briehy develop the general 
properties of such freeze-in scenarios, following Ref. |88j . 

Suppose that the coupling of FIMPs, X, to Standard-Model particles is 
characterized by a dimensionless number, A ;< 1, and the mass of the heaviest 
particle is denoted by m. Then the FIMPs are produced by interactions 
of the thermal bath of Standard-Model particles with a characteristic cross 
section that is on the order ax ~ A^/T^ for T > m. Since the density of 
Standard-Model particles is ~ T^, during one Hubble time, fn — Tfpi/T^, 
the physical density of FIMPs increases by Anx axT^tfj ~ A^AfpiT^. 
Number densities, n, are often expressed in terms of the yield, Y = n/s, 
where s ~ is the entropy density which is itself on the order of the photon 
number density. Yields are only changed by interactions and are unaffected 
by the expansion of the universe. In our case we thus have AYx ~ A^Mpi/T, 
which is largest at low temperatures and applies as long as the abundances 
are not Boltzmann-suppressed, and thus for T > m. One thus hnds 




m 


(28) 
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for the yield today. For comparison, the standard yield of thermal WIMPs 
is Yx — Pc^x/i'fnxSo), where pc and sq are the critical density and entropy 
density today, respectively. Using Eq. Q this gives 


- 

Mpimx(o-xx'^’} 


1 mx 

A2 Mpi ’ 


(29) 


where in the last expression we have phrased the cross section as {cxxv) ~ 
Identifying m with mx, one now realizes that Eqs. (28) and (29) 
are just the inverse of each other ! Thus, whereas the WIMP abundance 
decreases with increasing coupling. A, the FIMP abundance increases with 
A. Candidates for FIMPs include moduli or their supersymmetric partners 
in string theories, Dirac neutrinos within weak-scale supersymmetry, massive 
gauge bosons of an extra U{1) group that mix with photons, as already 
encountered in Sect. 2.2.2, and weakly coupled Kaluza-Klein states from 


extra dimensions as discussed in Sect. | 

Another alternative discussed in the literature is known as asymmetric 
dark matter. Usually it is assumed that the initial abundance of dark-matter 
particles, X, is equal to that of their anti-particles, X. One could, however, 
also imagine a situation in which X and X are distinguishable and initially 
not of equal abundance. In this case, in thermal freeze-out, for example, the 
more numerous component would always survive so that the relic abundance 
(cf. Eq. would only describe the symmetric component. The annihilation 
cross section could then be much larger, potentially permitting complete 
annihilation of the symmetric partner. In this situation the hnal dark-matter 
abundance would be determined by the initial asymmetry, analogous to the 
baryon number in certain scenarios of baryo- and leptogenesis. It has been 
speculated that unknown physics beyond the Standard Model could render 
comparable the hnal number densities for dark matter and baryons, so that 
flx/X.b — mx/mx, suggesting a dark-matter mass mx of a few GeV. If 
today dark matter would be constituted only by the symmetric part there 
would essentially be no remaining annihilations today, and thus the indirect 
signatures due to dark-matter annihilation that we will discuss in Sect. 
would be moot. The remaining dark matter could, however, still be trapped 
in astrophysical bodies on account of scattering with ordinary matter. 
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3. Direct searches 

3.1. Detection principles 

Direct-detection experiments rely on the scattering of dark-matter par¬ 
ticles from the halo of the Milky Way in a detector on Earth. The latter 
is usually set up deep underground, e.g. at Gran Sasso in Italy, Sanford in 
South Dakota or in the future Jin-Ping in China, to shield it from cosmic 
radiation, which, together with natural radioactivity in the rock and the de¬ 
tector material, constitutes the most important background. To discriminate 
the dark-matter signal from background processes and thermal noise, one ex¬ 
ploits usually two out of three possible signals: energy deposition (phonons) 
in (often cryogenic) calorimeters, scintillation light (photons), and ionisation 
(electrons). Localization of the scattering event allows one to define a fiducial 
volume inside the detector, which is then effectively shielded by the detector 
material itself. For the latter, current experiments use CaW 04 (CRESST 
HSHnniini]), CsI(Ti) (KIMS 1121193 ]), Ce (CDMS [HllSl [96l EHEHlIH] , Co- 
GeNT unnilini], edelweiss [Ml [103]), Nal(Ti) (DAMA/LIBRA [TnTj i 
or Si (CDMS Si [IT]) crystals, but also liquid Ar (DarkSide |105j i. Ne (DEAP 
[HE]) or Xe (LUX [TOT], XENON [MlIM 11101HII] , XMASS [^, ZEPLIN 
|113] ). For many of these experiments, upgrades are in preparation (e.g. 
DARWIN HU], EURECA [Il5], LZ [Il6], PANDA-X [HT], XENONlT/nT), 
in some cases focussing on the conversion of axions into photons (ADMX 
HH). It is also possible to exploit the metastability of a superheated liq¬ 
uid composed of carbon and fluorine (COUPP jll9j . PICASSO |12Uj . PICO 
HZH, SIMPLE |122j ). the annual modulation of the dark-matter signal as 
the Earth orbits around the Sun (DAMA/LIBRA [104] i. or the direction of 
the dark matter “wind”, blowing from the Cygnus constellation towards the 
Sun as it orbits around the galactic center (DMTPC |123] . DRIFT [124] . 
MIMAC HH], NEWAGE [HE] i. 

3.2. Local dark matter density and velocity 

The direct-detection rate depends on the local dark-matter density, cur¬ 
rently estimated to be p ~ 0.39 ± 0.03 GeV cm’ -3 mi- The systematic un¬ 
certainty of this value is up to 40%, on account of, e.g., effects of the galactic 
disk, whereas small halo substructures (clumps or streams) are unlikely to 
influence it signihcantly. It is therefore commonly adopted as a Standard 
Halo Model (SHM) value in halo-dependent comparisons of direct-detection 
experiments [I2H|. 
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The direct-detection rate also depends on the local dark-matter velocity 
distribution /(v). If we assume the velocity distribution at the solar circle to 
follow a Gaussian and take into account the orbit velocity of the Sun around 
the galactic center, vq = 220 km s“^, and the orbit speed of the Earth around 
the Sun, v^, we obtain the normalized one-dimensional Maxwell-Boltzmann 
distribution 


f{v)dv 


vdv 

VeVo^/T: 





— exp 


{v + 


(30) 


which must in addition be truncated at the local escape velocity from the 
Milky Way km s ^ ra- The dark-matter flux passing through 

a detector on Earth can therefore be quite large, nv = pv/irtx — 10^ cm“^ 
s~^ for mx = 100 GeV. 


3.3. Recoil energy and rate 

The recoil energy induced by a dark matter particle of mass mx on a 
nucleus of mass mx is 


Er 


2m X 


-(1 — cos 9*), 

mx 


(31) 


where q is the momentum transfer, /i = mxmx/{mx + mx) the reduced 
mass, V is now the mean dark-matter velocity relative to the target, and 9* 
the scattering angle in the center-of-mass system m- If we allow for inelastic 
scattering, i.e. that the dark-matter particle leaves the detector in an excited 
state, the minimal and maximal recoil energies 
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acquire a dependency on the increment S of the dark-matter rest mass, mx 
USD]. In the elastic case {S = 0) and for mx < vrix, fhe maximum 


^max 
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(33) 


ranges typically between 1 and (9(100) keV, the upper limit imposed by 
Uesc- Dark matter with masses below 1 GeV is typically below the detection 
threshold for nuclear recoil, but it could deposit enough energy to interact 



















with electrons (1 to 10 eV). It could then be detected via electron ionization, 
excitation or molecular dissociation. This option is particularly interesting 
for candidates preferably interacting with leptons, so-called leptophilic dark 
matter, that have received attention in the context of the observed positron 
excess (see below). 

In a detector composed of different target nuclei, the total recoil energy 
spectrum 
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is obtained by summing the differential rates for all target nuclei N, weighted 
by their number per unit mass of detector {Nx)- The individual nuclear recoil 
rates depend on the number density of the dark-matter particle, p/mx, which 
must be rescaled if the considered candidate is not the only contributor to 
the total relic density, and the differential scattering cross section 

^ MEr) (35) 

dEn 2p?v^v‘^ — 25/p, 

integrated over the velocity distribution with a lower limit of integration. 
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that depends on the minimal detectable recoil energy. The recoil rates are 
usually small, of order 1 event/100 kg/day. 

Experimentally, the recoil energy spectrum must still be convolved with 
an energy-resolution function and multiplied with an energy-dependent count¬ 
ing efficiency or cut acceptance. The probability that only a fraction of 
the true recoil energy is experimentally visible is incorporated in a so-called 
“Q-factor”, Q{Er) = E^is/Er. E^is is often quoted in keVee (keV electron- 
equivalent) or in photoelectrons. Since both background and signal events 
mostly involve scattering off electrons and nuclei, respectively, the Q-factor 
is separately determined from y-ray (e.g. ^'^Co or ^^^Cs) and neutron (e.g. 
AmBe or ^^^Cf) calibration sources. In crystals, where mostly phonons are 
produced, it ranges between 0.09 for iodine and 0.3 for germanium. In noble 
gases, a similar factor measures the scintillation efficiency of dark matter 
particles relative to photons. Both correction factors are subject to large 
experimental uncertainties at low energies. 
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3.4- Spin-independent cross sections 
Spin-independent (SI) cross sections 

<(£„) = lZ + (A-Z)(f„/f,)f^a,Fl,(Ea), ( 37 ) 

resnlt from scalar or vector conplings. Here, Z and A are the proton and 
total nncleon nnmbers of the target nnclens N, fp(^n}, hp ^ind cXp are the dark- 
matter conplings, rednced mass and total cross section for protons (nentrons), 
and Fsi{Epi) is the spin-independent nnclear form factor, i.e. the Fonrier 
transform of the nncleon density in the nnclens. It is often assnmed to be 
similar for protons and nentrons and have the Helm form 

Fi^iEn) = (38) 


with q = yj2mjqEpi, the effective nuclear radius tq ~ ■>/( 1.2 fm — 5 ^ 2 ^ 

and the nuclear skin thickness s ~ 1 fm m- For isospin-conserving cou¬ 
plings {fn = fp), the nuclear cross section is enhanced by a factor over 
the proton cross section, i.e. the whole nucleus interacts coherently provided 
the inverse momentum transfer is smaller than the size of the nucleus. Since 
heavier nuclei consist of more neutrons than protons, isospin-violating cou¬ 
pling with fn/fp — —Z/{A — Z) can reduce the cross sections for particular 
nuclei, and thus weaken the tension between conflicting experimental results 
|132] . To adjust the sensitivity of Xe- and Ge-based experiments to the signal 
regions claimed by CDMS-H, CoGeNT, and most strikingly DAMA/LIBRA 
(see below), one needs, e.g., fn/fp = —0.7 and —0.8, respectively. 

If, e.g., dark matter scatters off nuclei via the exchange of a scalar particle 
0 with mass and dimensionless coupling constants Ax, Ap and A^ to dark 
matter, protons and neutrons, respectively, the differential cross section is 

daf} _ ruN A^ [ZX^FpiEn) + (A - Z)XnFn{ER)f 

dEn ( 5 '^ + ^^)^ 27rn (n^ — 2(5//i)^^^ 

Since — Eff'^ = 2/i^n (n^ — 2(5//i)^^^/mx, the total spin-independent 

scattering cross section is then of the order 


a 


{Er) ~ 


iflniRrER + m2) 




{Er). 


(40) 
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Currently, three direct-detection experiments claim evidence for poten¬ 
tial signals of dark matter. At Gran Sasso, DAMA/LIBRA with 14 annual 
cycles and a total exposure of 1.33 ton-years hnds with 9.3cr conhdence an 
annual modulation of single-hit events in the (2-6) keV energy interval with 
amplitude (0.0112 ± 0.0012) counts/kg/keV/day, a measured phase of (144 
± 7) days and a period of (0.998 ± 0.002) years, all well in agreement with 
those expected for dark-matter particles |104j . CoGeNT at the Soudan Un¬ 
derground Laboratory claims a 2a (now somewhat less) dark-matter excess 
and annual modulation HU, and GDMS-II at SNOLAB observes three un¬ 
explained low-energy events in their Si-detector data sample. After their 
detector upgrade, GRESST-II no longer hnds a previously claimed excess 
[89] [91]. These claims are shown in Fig. (left) as colored areas. They 
cluster in the mass region of tens of GeV and at cross sections between 10“^^ 
and 10“^® cm^. All other direct detection searches have set exclusion limits 
on the SI dark matter-nucleus cross section that contradict the claims listed 
above. They are shown in Fig. (left) as full lines. The currently best ex¬ 
clusion bounds come from XENONIOO (soon to be updated by XENONIT) 
mi and LUX [107] . The signal claims also seem to be in conhict with the 
hrst PANDA-X results mg. 

Asymmetric dark matter models (green ovals) have been constructed with 
a view to resolve the apparent contradiction between claimed signals and ex¬ 
clusion limits, whereas null results from the LHG have driven supersymmetric 
(big red circle) and extra-dimensional models (blue oval) to high masses and 
low cross sections (see below). DARWIN is designed to probe the entire pa¬ 
rameter region for dark-matter masses above ~6 GeV, limited by the neutrino 
background (yellow region) [133] . Experiments based on the mK cryogenic 
technique, such as SuperGDMS [91] and EUREGA [115] . have access to lower 
WIMP masses. The low-mass region is emphasized in Fig. (right), where 
the claims, including the area favored by GRESST-II [89] (light blue), are 
shown together the exclusion curves. We include that from the GRESST-II 
low-threshold analysis of a single upgraded detector module [H] (red), that 
did not conhrm the previous excess. 

3.5. Spin-dependent cross sections 

Spin-dependent (SD) cross sections 

a%^{En) = 32fi^Gl[{Jp, + l)/JN][{Sp)ap + {S^)ar,]^Fi^{ER) (41) 
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Figure 6: Left: Spin-independent WIMP-nucleon scattering results |134L 1133) . Claims 
from the crystal-based experiments CDMS Si |97j, CoGeNT |101j . CRESST [89] and 
DAM A/LIBRA |104j are shown together with exclusion limits from the noble-gas ex¬ 
periments ZEPLIN-III inSI, XENONIO [TOH], XENONIOO [1101, and LUX [TUT]. Also 
shown are the projected sensitivities of DarkSide-50 |105j . LUX [107) . DEAP3600 [106] . 
XENONIT, DarkSide G2, XENONnT (similar sensitivity as LZ [m]) and DARWIN 
[114j . In the yellow shaded region coherent neutrino scattering becomes an important 
background. Right: Spin-independent WIMP-nucleon scattering results at low mass m- 
The 90% G.L. upper limit from GRESST-II (solid red) is depicted together with the 
expected sensitivity (Itr G.L.) from the background-only model (light red band). The 
GRESST-II 2a contour reported for phase 1 [SO] is shown in light blue. The dash-dotted 
red line refers to the reanalyzed data from the GRESST-II commissioning run [90]. Shown 
in green are the limits (90% G.L.) from Ge-based experiments SuperGDMS (solid) [00] . 
GDMSlite (dashed) [05] and EDELWEISS (dash-dotted) [102] . The parameter space fa¬ 
vored by GDMS Si [07] is shown in light green (90% G.L.), the one favored by GoGeNT 
(99% G.L. [M]) and DAMA/LIBRA |104) in yellow and orange. The exclusion curves 
from liquid Xe experiments (90% G.L.) are drawn in blue (solid for LUX [107] . dashed for 
XENONIOO jllOj l. Marked in grey is the limit for a background-free CaW04 experiment 
arising from coherent neutrino scattering, dominantly from solar neutrinos |I35j . 


are due to axial vector coupling. They depend on the nuclear spin the 
dark-matter coupling to protons and neutrons, Up^n, the nuclear form factor 
FsoiEn) and in particular on the expectation values {Sp^n) of the proton and 
neutron spin in the target nucleus [136111371113811139) . Since these are of 
order unity, SD cross sections are smaller than SI cross sections by a factor 
Alf, and the bounds on the latter are thus considerably better than those on 
the former. 

To compare results obtained with different target materials, a common 
practice is to report the cross section for the interaction with a single proton 
or neutron, assuming no coupling to the other type of nucleon. Calculations 
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Figure 7: Spin-dependent WIMP-nucleon scattering results [lllj . XENONIOO 90%C.L. 
upper limits on the cross sections on protons (left) and neutrons (right) based on Ref. |139j 
(dark blue curve) are shown together with their la and 2a uncertainties on the expected 
sensitivity (green and yellow shaded bands) and results from XENONIO |109j based on 
Ref. dMl, CDMS [MlEH, ZEPLIN-III [m] based on Ref. [HH], PICASSO [IM], COUPP 
SIMPLE |122j . KIMS |M], and results derived from indirect IceCube limits |140j in 
the hard {W'^W~, and soft {hh) annihilation channels. 


of the structure functions, Sp^n, are traditionally based on the nuclear-shell 
model [1361 I137|, 1138] . but have recently also been performed using chiral 
effective held theory, yielding much better agreement between calculated and 
measured Xe spectra, both in terms of energy and in ordering of the nuclear 
levels |139j . The currently best upper limits from XENONIOO for protons and 
neutrons are shown as full blue curves in Fig. (left and right) with their la 
and 2a uncertainty bands (green and yellow), together with results from other 
competing experiments. The sensitivity to proton couplings is weaker, since 
both Xe isotopes (^^®Xe and ^^^Xe) have an unpaired neutron, but an even 
number of protons leading to |(5'p)| |(5'„)|. The IceCube detector at the 

south pole has searched for muon neutrinos from dark-matter annihilation 
at the center of the Sun. This indirect-detection mode will be discussed 
in Sect. |4.5[ If one assumes a particular dark-matter candidate, such as 


the lightest neutralino of supersymmetry |141j . the indirect-detection limits 
from heavy quark (6), lepton (r) and weak gauge boson iW^) annihilation 
products can be translated into competitive SD direct-detection limits [140] . 
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3.6. Light dark matter 

As discussed in Sect. axions were introduced by Peccei and Quinn as a 
solution of the strong CP problem in the form of a pseudo-Nambu-Goldstone 
boson emerging from the breaking of a global U(l) symmetry [85]. Although 
the original model has been ruled out, invisible axions arising from a higher 
symmetry-breaking scale are still allowed, as explicitly demonstrated in the 
DFSZ and KSVZ models [HSl EM EM EM- 

The Sun is believed to be an intense source of axions through produc¬ 
tion via bremsstrahlung, Compton scattering, axio-recombination and axio- 
deexcitation which are standard electrodynamic processes, in which the hnal- 
state photon is replaced by an axion |146j . Conversely, axion-like particles 
(ALPs) could have been generated via a non-thermal production mechanism 
in the early universe, so that they would now constitute the dark matter in 
our galaxy. Both types of particles give rise to direct observable signatures 
through their coupling to photons (qa'i), electrons {gAe), and nuclei {gAN)- 

The coupling gAe niay be tested via scattering of electrons through the 
axio-electric effect [1471114811149[ 115011151j . Fig. (1 eft) shows the current 
exclusion limits on g^g. The XENONIOO results (90% C.L., blue line) are 
shown together with the expected sensitivity at Icr and 2a based on the 
background hypothesis (green and yellow bands). The range of validity of the 
light solar-axion flux limits the search to rriA < 1 keV |146j . For comparison, 
other recent direct search limits and astrophysical bounds are also shown, 
together with the theoretical benchmark models DFSZ and KSVZ. For solar 
axions with masses below 1 keV, the electron coupling can be limited by 
the XENONIOO data to gAe < 7.7 x 10“^^ at 90% C.L. Within the DFSZ 
and KSVZ models, this limit excludes axion masses above 0.3 eV and 80 eV, 
respectively. For comparison, the CAST experiment has tested the coupling 
to photons QA-y and excluded axions within the KSVZ model in the mass 
range between 0.64 eV and 1.17 eV [TS2]- 

The current exclusion limits on galactic ALPs from XENONIOO, compet¬ 
ing experiments and astrophysical bounds are shown in Fig. (right) together 
with the KSVZ benchmark model. The steps in the sensitivity around 5 and 
35 keV reflect jumps in the photoelectric cross section reflecting atomic en¬ 
ergy levels. In the 5-10 keV mass range, XENONIOO sets the best upper 
limit, excluding an axion-electron coupling gAe > 1 X 10“^^ at the 90% C.L., 
if one assumes that ALPs constitute all of the galactic dark matter. This 
case has been discussed in Sect. 
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Figure 8: Left: Current limits on solar axions |153j . The XENONIOO limits (90% C.L., 
blue line) are shown together with the expected sensitivity at la and 2a based on the 
background hypothesis (green and yellow bands). Also shown are limits by EDELWEISS- 
II |103j and XMASS |II2) together with those from a Si(Li) detector |I54) . Indirect 
astrophysical bounds from solar neutrinos |155j and red giants |I56) are represented by 
dashed lines. The benchmark DFSZ and KSVZ models are represented by black solid 
lines |142111431II441IT^ . Right: Current limits on galactic axions |153j . The XENONIOO 
limit (90% C.L.) on ALPs coupling to electrons as a function of the mass is shown under 
the assumption that ALPs constitute all the dark matter in our galaxy (blue line). The 
expected sensitivity is shown by the green/yellow bands (la/2a). The other curves are 
constraints set by CoGeNT |10()j (brown dashed line), CDMS |HS] (grey continuous line), 
and EDELWEISS-II |103] (red line, extending up to 40 keV). Indirect astrophysical bounds 
from solar neutrinos |155j are represented as a dashed line. The benchmark KSVZ model 
is represented by a black solid line |144[ 1145) . 

4. Indirect searches 

4.1. General facts on indirect detection of WIMPs 

Dark matter cannot only be detected directly in dedicated experiments 
searching for nnclear recoils from the scattering of dark-matter particles or 
prodnced in particle accelerators snch as the LHC, bnt it can also reveal 
its existence indirectly. The total nnmber of dark-matter particles does not 
change signihcantly after freeze-ont in the early nniverse, bnt their spatial 
distribution changes considerably during structure formation. The very self- 
annihilation that plays a central role in this freeze-out can give rise to a 
signihcant flux of y-rays, neutrinos, and even antimatter such as antiprotons 
and positrons, especially in regions with large dark-matter density. The 
energy of the secondary particles can reach up to the mass of the dark-matter 
particle, which typically would be a few hundred GeV. 
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Positrons would eventually annihilate with the electrons of the interstellar 
plasma and give rise to 511-keV line emission. In fact, such a line has already 
been observed with INTEGRAL |157j . in particular from the galactic center, 
but uncertainties remain as positrons are copiously produced by [3 decay of 
unstable isotopes Secondary electrons and positrons can give rise to 

synchrotron radiation that may be detected in the radio band. Therefore, 
detectors for cosmic rays and 7-rays, neutrino telescopes, and radio telescopes 
can be used for indirect dark-matter detection. 

Since dark-matter annihilation scales with the square of its density, indi¬ 
rect detection is more sensitive to cosmological and astrophysical processes 
than is direct detection. Such processes include the appearance of cusps 
around the central super-massive black holes in our own as well as other 
galaxies and the formation of clumps that condense around intermediate- 
mass black holes or that arise from initial density fluctuations. 

As the fluxes and spectra of ordinary cosmic rays are modihed during 
propagation to Earth, so are those of the annihilation and decay products 
of dark matter. Charged particles undergo diffusion in the galactic magnetic 
field, whereas photons and neutrinos will propagate on straight trajectories 
and can therefore reveal the spatial distribution of dark matter. Eventually, 
this might lead to a measurement of the small-scale structure of dark matter, 
thus permitting an experimental test of our theoretical understanding of the 
formation and survival of cusps and clumps. 

The annihilation yield of a dark-matter particle X is usually calculated 
under the assumption that the particle in question is of Majorana type, i.e. 
it is its own antiparticle, X. Otherwise there would be an additional factor 
1/2. Then, the flux of neutral secondaries can be written as 

E) = ^ dl A !.■((. n)l. (42) 

where the index i denotes the secondary particle observed (we focus on 7-rays 
and neutrinos) and the integration variable I refers to the path length along 
the line of sight (l.o.s.) in direction n. Furthermore, dNi/dE is the multiplic¬ 
ity spectrum of secondary particle i, mx is the dark matter particle mass, 
and n)] is the dark-matter density along the line of sight. One often 

assumes that the annihilation cross section averaged over the dark-matter 
velocity distribution, (cxxv), is the same as that relevant for thermal relic 
freeze-out, i.e. Eq. Note, however, that the particle velocities relevant for 
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indirect detection are very different from those pertaining to thermal freeze- 
out: the typical WIMP velocity in the Milky Way is n ~ 300 kms“^, whereas 
freeze-out happened when Tf ~ mx/20, implying v ~ 0.4 c. Therefore, if the 
WIMP annihilation cross section has a strong velocity dependence, as is 
possible if S'-waves are suppressed and higher partial waves dominate, then 
constraints from freeze-out and from indirect detection may not be directly 
comparable. 

To separate the factors depending on the dark-matter prohle from those 
depending only on particle physics, we introduce, following |159j . the quantity 

■'<"^ = 8l^( o.3GeVcm-3 ) L.s. 

where we used galactic scales and px — 0.3 GeVcm”^ to render J dimension¬ 
less. Note that various dehnitions of J can be found in the literature. We 
dehne J(n, AG) as the average of J(n) over a circular region of solid angle 
AG, centered around n. 

The parameter J describes the intensity distribution of neutral secon¬ 
daries, whereas J is a measure of the flux from the region AG, given by 


$i(AG,E) ~ C 


dNi 

dE 


pb J 


/lTeV\ 
V mx y 


J(AG) AG 


(44) 


with C = 2.75 X lO'^^ cm-^s-P 

The enhancement of annihilation rates due to small-scale structure in the 
form of dark-matter clumps (cf. Eq. is often expressed in terms of a so- 
called boosting factor, B ~ (p^)/(p)^, that measures the annihilation yield 
relative to that computed using a smooth dark-matter distribution. 

Another type of signal amplihcation can arise, if dark-matter interac¬ 
tions involve as a force carrier a new boson, 0, of sufficiently low mass, 
Mrf) <10 GeV. Originally envisioned by Sommerfeld in a different context 
preu] . the wave function of the incoming particle, conventionally approx¬ 
imated as a plane wave, can be signihcantly distorted, leading to strong 
upward or downward variations in the cross section. If the parameters are 
chosen appropriately, this so-called Sommerfeld enhancement can lead to a 
substantial increase in the annihilation cross section in galaxies and clusters, 
while leaving the arguments about thermal freeze-out in the early universe 
intact m- Furthermore, certain lowest-order cross sections can be sup¬ 
pressed by symmetries or helicity effects, such as the annihilation of two 
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non-relativistic Majorana fermions into a relativistic fermion/anti-fermion 
pair. In this case internal bremsstrahlung^ the attachment of an additional 
photon to one of the charged hnal-state fermions, can actually enhance the 
annihilation cross section, because the larger available phase space alleviates 
suppression arising from symmetry or helicity effects. 

If, instead of annihilation, dark matter dominantly decays with lifetime tx 
and a multiplicity spectrum of secondary particles dNi/dE, the flux observed 
on Earth would read 

$i(n,E) = ^- f dl px [r(/,n)] . (45) 

dE ATiTxmx Jlo.s. 

Note that in this case there is no boosting, because the flux depends only 
linearly on the dark-matter density. 


4-2. Photons 

Let us begin with a possible signature of dark matter in the X-ray sky. Re¬ 
cent analyses of the X-ray spectrum observed toward the Andromeda galaxy, 
the Perseus galaxy cluster, and various other sources stacked together in¬ 
dicated a 3.55-keV line that is not expected as a thermal emission of an 
astrophysical plasma [HI 1162] . If interpreted in terms of decaying sterile 
neutrinos, one derives a neutrino mass ~ 7.1 keV and a mixing angle 
with active neutrinos sin^ 9g = (5.5 — 50) x 10“^^. In Fig. [^of Sect. these 
values are indicated by a red square and found not in conflict with other 
constraints. 

The detection of an X-ray line and its interpretation is the subject of an 
ongoing debate, though. The Suzaku satellite hnds a line in the X-ray emis¬ 
sion of the Perseus galaxy cluster, but not for other galaxy clusters usa. 
Concerns have been raised that the intensity distribution of line emission 
doesn’t correspond to that expected for dark-matter decay [TM] . No X- 
ray lines are observed from dwarf galaxies |165] . for which confusion with 
thermal line emission should be absent, and other galaxies [166] . Several 
theoretical pre- and postdictions are related to this possible X-ray line, in¬ 
volving axion-like particles UST], axion-photon conversions |168j , that we will 
discuss in general in Sect. 4.6, and scenarios in which inelastic scattering of 


dark matter leads to an excited state that subsequently decays |169] . Among 
alternative astrophysical explanations, the well-known atomic lines K XVIII 
or Cl XVII have been proposed |170] and criticized on account of large sys¬ 
tematic uncertainties jl71j or wrong transition probabilities in the case of Cl 

XVII [TT2] . 
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Figure 9: Spectral energy density, E'^ x AN/dE, of 7-rays from the galactic-center source 
seen by H.E.S.S. m, complemented with data points obtained with VERITAS [mi- 
upper limits are 95% CL. The shaded area shows a power-law fit, dN/dE ^ E 
The dashed line illustrates typical spectra of phenomenological minimal supersymmet¬ 
ric standard-model annihilation for neutralino masses of 14 TeV. The dotted line shows 
the distribution predicted for Kaluza-Klein dark matter with a mass of 5 TeV. The solid 
line gives the spectrum of a 10 TeV dark-matter particle annihilating into r+T“ (30%) 
and bb (70%). 


Let us now turn to 7 rays. Indirect detection of WIMPs requires disentan¬ 
gling the flux of dark-matter annihilation products from that of conventional 
astrophysical contributions. Several observations have been interpreted as 
possible dark-matter annihilation signatures. These include TeV-band 7 - 
ray emission imi from the galactic center, described in Fig. and the 
INTEGRAL observation of a 511-keV 7 -ray line from the galactic-center re¬ 
gion m, elucidated in Fig. 

However, with more and newer data, these signatures seem to find con¬ 
vincing explanations involving normal astrophysical processes: The TeV 7 - 
rays seen by H.E.S.S. and VERITAS |174j display no spectral structure up 
to > 30 TeV which would require an unnaturally heavy dark-matter primary. 
Instead, the spectrum appears compatible with predictions based on an ac¬ 
celerated primary cosmic-ray component |175] . In addition, even for the 
cuspy NFW profile, the required velocity-weighted annihilation cross section 
into 7 rays would be of the order 3 x 10“^^cm^s“^, higher by a factor 100 
than that needed to produce the correct thermal relic abundance, cf. Eq. ([^. 
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Figure 10: The 511-keV emission seen with INTEGRAL/SPI. Left panel: The spectrum 
from the inner galaxy (b < 15°) compared to fits to a combination of astrophysical com¬ 
ponents. Right panel: Longitudinal distribution of line flux compared to predictions for 
dark-matter annihilation for various density profiles (right panel). Taken from Ref. |158) . 


Whatever the interpretation, the measurement provides an upper limit to the 
annihilation cross section nisi. The 511-keV intensity distribution seen with 
INTEGRAL does not show the spherical symmetry expected for the central 
region of a dark-matter halo, but looks more like the galactic bulge nm. 




Figure 11: The Fermi GeV y-ray excess. Left panel: The residual y-ray spectrum after 
subtraction of astrophysical y-ray emission averaged over the region 2° < |6| < 20° and 
|1| < 20°. Also shown are best-fit model spectra of two dark-matter annihilation channels. 
From Ref. |182| . Right panel: Intensity profile of the excess at = 2 GeV compared with 
a contracted NFW profile oc From Ref. |18.8j . 

In recent years, the LAT detector aboard the Fermi satellit^has surveyed 
the 7 -ray emission from the inner galaxy. One result is evidence of excess 
emission extending up to 40° above and below the galactic plane, that has a 


^http://www-glast.Stanford.edu/ and http://fermi.gsfc.nasa.gov/ 
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spectrum roughly following an E~‘^ power law, significantly harder than that 
expected of emission from the galactic disk |178] . This structure is collo¬ 
quially known as “Fermi bubbles” and roughly matches a radio counterpart 
found earlier, the so-called “WMAP haze” |179] that was recently conhrmed 
by the Planck experiment HHDl. The WMAP haze has a smaller vertical ex¬ 
tent than the Fermi bubbles, possibly on account of the difference in the ra¬ 
diative energy-loss time of synchrotron-radiating and inverse-Compton scat¬ 
tering electrons. Whereas the Fermi bubbles may be interpreted in terms of 
outflows triggered by past black-hole activity or a central starburst, a second 
7 -ray excess on smaller scales has been claimed that may be related to dark 
matter |181j . In Fig. 11 we show the spectrum and latitudinal distribution 
of this GeV excess. If it is interpreted in terms of dark matter, one would 
require a particle mass in the range of 7 — 12 GeV for leptonic channels or 
25 — 45 GeV for hadronic hnal states. The cross section would be roughly 
consistent |184] with that estimated from thermal freeze-out, cf. Eq. Q, 
although its value would be uncertain by up to a factor 50 on account of 
the poorly known density profile of dark matter. In any case, conventional 
astrophysical scenarios cannot be excluded, and so one can, strictly speaking, 
only place an upper limit on the cross section of dark-matter annihilation. 

A word on y-ray producing annihilation channels is in order. Depending 
on the nature of the initial annihilation products, two processes are to be din- 
stinguished. First, a pair of quarks may be created either directly or through 
electroweak gauge boson pairs. The quarks subsequently evolve into, e.g., 
neutral pions that decay into 7 rays. Second, charged annihilation products 
should emit bremsstrahlung. Both lead to a broad energy distribution of 7 
rays that is easy to confuse with contributions from other astrophysical pro¬ 
cesses. Truly incontrovertible evidence for dark matter would be a 7 -ray line 
because we do not know an astrophysical process that would produce a line 
at energies > 1 GeV. Unfortunately, the process XV —)■ 77 is not possible 
with tree-level processes because dark matter is neutral, and it requires a 
loop of virtual charged particles. 

As GDM WIMPs are by definition slow, annihilation into two photons 
will inevitably give a line at E.y = mx, and so one would automatically mea¬ 
sure the particle mass as well. The channel XX —)■ 7 Z, H would also lead 
to a narrow 7 -ray line, albeit at the energy E^ = mx [l — nil „/(4m^)]. 
Recent reports of a 7 -ray line around 130 GeV |185[ 1186] caused substan¬ 
tial excitement, but subsequent analysis suggested that the feature is likely 
caused by either a statistical or an instrumental effect |187] . 
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Among all targets of 7 -ray observations, dwarf spheroidal galaxies are 
of particular interest, because they are dark-matter dominated and thus 
less subject to uncertainties from astrophysical flux contributions. Upper 
limits have been derived by both ground- and space-based 7 -ray detectors 
The coincidence of a 7 -ray source with the dwarf galaxy 
Reticulum 2 is now subject of further investigation |192j . Other constraints 
have been derived from the galactic halo and subhaloes, galaxy clusters (for 
a recent summary see, e.g., |193] ). and from the extragalactic 7 -ray back¬ 
ground |194j . For mx ^ 100 GeV, 7 -ray observations challenge the thermal 
freeze-out scenario, cf. Eq. and the thermal-relic model is excluded for par¬ 
ticle masses below mx ~ 10 GeV. The upcoming GTA experiment|^ promises 
to provide the sensitivity to test the freeze-out model up to the TeV-scale. 


4-3. Electromagnetic cascades and their effects on the CMB 

Dark-matter annihilation with a non-vanishing branching ratio into the 
electromagnetic channel leads to distortions of the GMB which can be probed 
with data of, e.g., WMAP or Planck. The analysis is simplihed, if the annihi¬ 
lation products are relativistic. In that case, the rate of change in comoving 
energy density can be estimated as 


d [ p 

~dt [(TT^ 


mx {(Jxxv) 


nxi^) 

{l + zY 


c? 


mx 


{l + zf. 


(46) 


GMB observations are sensitive to the injection of additional energy, and 
so the absence of signatures of additional energy sources provides an upper 
limit on the annihilation cross section, for which the conversion efficiency 
of turning energy into GMB distortions is parametrized by a factor /efr(^). 
Results of such an analysis are summarized in Fig. 12 [21|. To be noted from 


the hgure is the exclusion of the thermal freeze-out scenario for low WIMP 
masses. The dependence of the upper limits on the WIMP mass is linear, as 
it is for the 7 -ray limits, and the advantage of GMB-based limits lies in the 
absence of astrophysical uncertainties. The efficiency parameter fesiz) has 
a relatively small uncertainty. 

If we assume that the thermal-relic scenario is essentially correct, then 
the combination of the corresponding cross section (cf. Eq. with the 


^http://www.eta-observatory.org/ 
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Figure 12: Constraints on the self-annihilation cross section at the time of recombination, 
{av)z,, times an efficiency parameter, /eg. The blue area shows the parameter space ex¬ 
cluded by the TT, TE, EE, and low-P data of Planck at 95 % CL. The yellow line indicates 
the constraint using WMAP9 data. The dashed green line delineates the region ultimately 
accessible by a cosmic-variance-limited experiment with angular resolution comparable to 
that of Planck. The horizontal red band denotes the thermal-relic cross section scaled 
with /eff for different annihilation channels. The dark grey circles show the best-fit DM 
models for the cosmic-ray positron excess |195) . The light grey stars refer to best-ht mod¬ 
els for the galactic-center gamma-ray excess m. with the light grey area indicating the 
astrophysical uncertainties. Taken from [M| . 


constraints displayed in Fig. 12 provide a lower limit on the WIMP mass 

[e.g. me], 

mx > mcMB = (100 GeV) fes{Zrec) ■ (47) 

On the other hand, the same assnmption combined with cross section nni- 
tarity, (cr^xf) ^ l/m^, implies mx < 20TeV. One would then conclude 
that thermal-relic WIMPs with signihcant branching to electromagnetic an¬ 
nihilation products must have a mass between a few tens of GeV and about 
20 TeV. 

If the WIMP (X) carries only part of the dark matter density, Ox < 
then the constraint must be adjusted. If we denote with (cxth'^) the thermal- 
relic cross section saturating and (o'cmb'?^) is the GMB constraint (the 
thick blue line in Fig. 12), then the cross section limits change to {(Txxv) > 
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{athv) for mx > mcMB and {axxv) > (cTth^^) V(<^cmb^^) for mx < rricuB- 
4■4- Antimatter as a product of dark-matter interactions 




Figure 13: Left panel: Compilation of data on the GeV “excess” of galactic positrons. 
Expectations for interstellar production and speculations on, e.g., a contribution from 
pulsar-wind nebulae are added. Right panel: The electron plus positron flux as measured 
with various detectors. The shaded bands indicate the uncertainties arising from modeling 
contributions of primary electrons and positrons from supernova remnants, pulsar wind 
nebulae, and from secondary electrons and positrons produced during propagation. Taken 
from Di Mauro et al. m- 


There is now strong evidence for an increase with energy of the cosmic-ray 
positron fraction, + N ~). In Fig. 13 we summarize results obtained 


with the PAMELA satellite |197] . Fermi LAT |198j . and AMS-02 |199j . that 
together clearly demonstrate an excess between a few GeV and a few hundred 
GeV, compared to the inevitable contribution from cosmic-ray interactions 
in the interstellar medium. This strengthened the case for a positron fraction 
in the electron plus positron flux that is increasing between ~ 10 GeV and 
~ 100 GeV (see Fig. 13, left panel). Gonventional wisdom is that all positrons 
are secondary particles produced through collisions of cosmic rays with gas in 
the interstellar medium. As a function of energy, the interstellar contribution 
to the positron fraction should fall off. The observed positron excess there¬ 
fore calls for additional sources of positrons, a candidate for which would be 
pulsars or pulsar-wind nebulae. Both primary and secondary positrons can 
be produced within the sources by electromagnetic cascades induced by accel¬ 
erated electrons [200]. Interactions of accelerated primary particles |201] are 
possible, but imply similar structures in the spectra of secondary cosmic-ray 
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elements |2U2j . In any case, fits of such astrophysical models can reproduce 
the rising positron fraction. 

PAMELA has not detected any signihcant enhancement of the antipro¬ 
ton flux beyond what is expected from interactions of cosmic rays with the 
gas in our galaxy |2U3] . Any interpretation of the positron excess in terms 
of annihilating or decaying dark matter thus requires so-called leptophilic 
dark matter, i.e. predominantly leptons as end products. In addition, the 
thermal-relic cross section (cf. Eq. is too small by a factor > 100 to pro¬ 
vide the observed positron flux. Correspondingly strong boosting would be 
required by either dark-matter clumping or Sommerfeld enhancement. In an 
interpretation involving dark-matter decay, the limits on the particle lifetime 
would be weak. 

The newest results by AMS-02 |2041I205] extend up to 500 GeV in positron 
energy and show that the positron fraction levels off and remains approxi¬ 
mately constant at energies above 200 GeV. If the positron excess is indeed 
due to dark matter, a particle mass of a few hundred GeV would be most 
likely. The implied extragalactic photon emission from electrons and positron 
severely challenge both scenarios, dark-matter decay |56] and annihilation 
|206j . Galactic photon emission constrains the excess production of electrons 
and positrons as well. Besides y-ray emission, synchrotron radiation in the 
radio band provides constraints on the branching of dark-matter annihilation 
into electrons and positrons. In particular for the galactic-center region, they 
are comparable to upper limits based on y-ray emission or possibly stronger, 
depending on the scenario [207] . Gare must be exercised, though, because 
the synchrotron emissivity does not only depend on the uncertain density 
prohle of dark matter in the innermost regions of the galaxy, as does the 
y-ray emissivity, it also hinges on the poorly known magnetic-held strength 
close to the galactic center. 

In principle, individual anti-nucleons produced by dark-matter annihila¬ 
tion or decay can coalesce to heavier anti-nuclei. Very little astrophysical 
background is expected for mass number A > 3. In addition, the back¬ 
ground from cosmic-ray interactions should also quickly fall off with increas¬ 
ing energy, whereas, e.g., the anti-deuteron flux from cold dark matter has 
a flat spectrum. No signihcant limits on the huxes of anti-nuclei have been 
placed to date. The planned General Antiparticle Spectrometer experimenlj^ 


“http://gammaO.astro.ucla.edu/gaps/ 
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(GAPS) |2U8j promises to provide constraints at kinetic energies np to 0.3 
GeV per nncleon, covering anti-protons as well, and the energy range aronnd 
1 GeV per nncleon shonld be accessible with AMS-02. 

In the right panel of Fig. [T^ we have already displayed the total electron- 
plns-positron spectrnm observed with Fermi LAT |2U9] . Measnrements with 
H.E.S.S. indicate a cnt-off or spectral break at abont ITeV |21Uj . Between 
~ 100 GeV and ~ 1000 GeV, an excess is visible compared to that expected 
for primary electrons from sonrces with nniform power-law spectra and ho- 
mogeneons spatial distribntion. Both assnmptions are not well jnstified. In 
fact, the observed spectrnm can be well explained by combining contribn- 
tions from individnal snpernova remnants and pnlsars. Fntnre experiments, 
that will measnre the combined electron and positron flnx np to several 
TeV, inclnde the balloon-based Gosmic Ray Electron Synchrotron Telescope 
(GREST) |212] and the GALorimetric Electron Telescope (GALET) |213j on 
the International Space Station. 

Model predictions of the electron and positron flnxes are obtained by solv¬ 
ing a diffnsion-energy loss eqnation that can be derived from a Fokker-Planck 
description of particle transport. Among the input parameters are the dif- 
fnsion coefficient and the relevant distribntions of both astrophysical sonrces 
and sites of dark-matter annihilation or decay. Nnmerical solntions can be 
obtained with dedicated software packages snch as GalPropj^and DRAGOhQ 
whereas analytical methods typically involve integrating the Green’s fnnction 
for the problem. 

The anisotropy of the galactic electrons and positrons can also serve as a 
probe for the origin of the excess positrons. The positron sonrce rate from 
dark matter is expected to be rather smoothly distribnted, unless strong 
boosting is dominated by high-mass dark-matter clumps, and thus a rela¬ 
tively small anisotropy is expected |214] . More detailed investigations predict 
an upper limit on the dipolar anisotropy given by 

l<5(£.)l < 3 X 10-=> for 30GeV < E, < 2 X lO’GeV. (48) 

This is consistent with current upper limits which are at the level of a few 
percent in this energy range |215j . The electron anisotropy for conventional 


® http: //galprop. stELnf ord. edu/ 

'http://www.desy.de/-macclone/DRAGON/index.html 
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sources such as SNRs is highly dependent on the actual location and activity 
period of the sources, but on average larger than that quoted in Eq. ^ by 
a factor of a few izig. Future detectors will have sufficient sensitivity to 
detect anisotropies of that order. 

Distinguishing dark-matter contributions from ordinary astrophysical sig¬ 
natures in the spectra of cosmic rays, antimatter, radiation fields from the 
radio to the 7 -ray band, and also neutrinos requires a multi-messenger ap¬ 
proach that includes accurate modeling of the acceleration and propagation 
of cosmic rays, as well as their production of secondary 7 -rays and neutrinos. 


Bergstrom, Bringmann & Edsjo (2010) 
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Figure 14: Experimental coverage of the cross-section plane of spin-independent WIMP- 
nucleon scattering and WIMP self-annihilation. As tracer of the latter we plot on the 
cc—axis {o'xx'^) /^x (‘-f- ^'l- 42) and restrict to 7 -ray observations of the galactic center. 
The density of dark matter profile is assumed to follow an NFW profile without additional 
boosting. All cross section combinations permitted in typical MSSM and mSUGRA models 
are shown. DMA denotes a possible dedicated dark matter array. From Ref. 


To be noted is the complementarity of direct and indirect dark-matter 


searches. In Fig. 14 we indicate the reach of current and future experiments 
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in constraining MSSM and mSUGRA dark-matter models through measure¬ 
ments of spin-independent WIMP-nucleon scattering and WIMP annihila¬ 
tion, the latter on the example of 7 -rays from the galactic center. The top 
half of the plot refers to large scattering cross sections addressed by the direct- 
detection experiments, whereas indirect detection probes the right half of the 
plot, where annihilation cross sections are large. It is interesting to note that 
Fermi LAT in combination with CTA will probe WIMP annihilation into 7 
rays down to very small values of the annihilation cross section |218j . which 
includes the region of small scattering cross section, and so CTA may reach 
deeper than the largest future direct-detection experiments. The mass range 
nix ^ 10 GeV will be difficult to probe both directly, on account of the small 
nuclear recoil, and indirectly, because low-energy secondary 7 -rays and neu¬ 
trinos would fall below the energy threshold of ground-based detectors. Such 
low-mass WIMPs are easier to probe at accelerator experiments, as we shall 
discuss in Sect. [H 


4-5. WIMP capture and neutrinos 

WIMP dark matter can be trapped in the gravitational well of celestial 
bodies such as stars and planets, resulting in an enhanced dark-matter den¬ 
sity near the center of these objects. An inevitable consequence is a large 
annihilation rate that depends on the WIMP-proton scattering cross section, 
as the following considerations illustrate. 

The rate of change of the number of WIMPs, Nx, is given by 

Nx = Tc-Ta = Pc , (49) 

Vx 


where Pc is the WIMP capture rate. Pa = CaN^ is the WIMP annihilation 
rate in the compact body, and Vx is the effective volume in which most 
of the WIMPs reside, typically the central core of the object in question. 
The steady-state solution to Eq. 49 and the relaxation timescale toward the 
equilibrium situation are 


A^eq 


TcEx 






X 


Tc {<^xx'^) 


(50) 


Obviously, knowledge of the capture rate. Pc, is required for estimating either 
quantity. A WIMP will be captured if it is scattered to a velocity smaller than 
the escape velocity of the astronomical body in question. Typical WIMP 
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velocities in the Milky Way, v ~ 300kms“^, are smaller than the escape 
velocity of the Sun, Vq ~ 620kms“^, and so capturing WIMPs in the Sun 
is an efficient process, that is dominated by axial (spin-dependent) WIMP 
scattering on hydrogen, the most common element in the Sun. The rate can 
be estimated as [T3] 


c,0 


~ (10 




Px 


L0.3GeVcm-3J 


1 

'100 GeV 


\ 1 
p 


'270 kms-3' 


mx 


10 “^° cm^ 


V 


(51) 

where <7®° is the spin-dependent WIMP-proton scattering cross section and 
the typical WIMP velocity v is comparable to the escape velocity of the 
galaxy. The escape velocity of the Earth, ^0 ~ llkms“^, is much smaller 
than the typical WIMP velocity, and so capture in the Earth has a very low 
efficiency. 


Using the parameter scales given in Eq. and further 
{axxv) ^ 10 




and Ux ~ 10 


27 cm3 


the relaxation timescale would be 


T., ss 10‘'s ~ 3- lO'yr. 


(52) 


For a wide range of WIMP masses and cross sections one would therefore 
expect a dark-matter density in the core of the Sun that is close to the 
steady-state value, whose ratio to the local dark-matter density, px/'^X) can 
be written as 


^eqmx ^ ^ q 12 

Vxpx 


(53) 


In the steady state, the total rate of annihilation events by dehnition equals 
the capture rate. Therefore, in this limit any secondary particle flux resulting 
from WIMP annihilation in the star depends on the same cross section that 
is probed by direct-detection experiments. In that respect, indirect detection 
by capture in the Sun is complementary to direct dark-matter detection. 

All hnal-state particles except neutrinos would be thermalized inside the 
star. The neutrinos would escape and can thus potentially be detected 
with large high-energy neutrino detectors such as IceCube. Of particular 
interest is the mass range mx > mw±, for which the annihilation channel 
X + X ^ W~^ + W~ is open and the subsequent decays provide a peak 
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in the neutrino spectrum near the energy mx/2. IceCube and its Deep- 
Core subarray have placed severe constraints on the cross section of both 
spin-independent and spin-dependent scattering m- Whereas for the spin- 
independent WIMP-nucleon cross section the IceCube limit at WIMP masses 
of a few hundred GeV and is not competitive with that of XENONIOO limits, 
the spin-dependent cross section is constrained to ;< 10“^° cm^ at WIMP 
masses of a few hundred GeV, which is the best limit published to date for 
mx ^ 35 GeV. It is also of the same order as cross sections predicted within 
the MSSM, so that the IceGube measurements permit to test that scenario. 

The Super-Kamiokande experiment has a much lower energy threshold 
compared to IceGube, which permits addressing the mass range 5 GeV < 
mx ^ 50 GeV. Preliminary upper limits on are on the order of 10“^® cm^ 
and 10“^° cm^ for annihilation into bb and rr, respectively. The future Hyper- 
Kamiokande and IceGube-PINGU experiments promise to reach sensitivities 
down to afP ~ cm^ between 10 GeV and a few hundred GeV. 


4-6. Astrophysical and cosmological signatures of WISPs 

In this section we consider signatures of WISPs in astrophysical and cos¬ 
mological observations that arise from the coupling of WISPs to photons. 
The physical concepts behind WISPs were described in Sect. |2.2[ The most 


well-known examples of WISPs are pseudo-scalar axions, ALPs, and hidden 
photons. In the presence of electromagnetic helds, in particular of magnetic 
helds, photons can oscillate into axions or ALPs and vice-versa, an effect 
known as the Primakoff effect IZE], 

If the wavelength of the photon and the ALP are much smaller than 
the length scales on which the properties of the medium and the magnetic 
held vary, the coupling of photons and ALPs is mathematically analogous 
to that of neutrino oscillations. The length scale of oscillations is typically 
much larger than those accessible with laboratory experiments, but orders of 
magnitude smaller than astrophysical and cosmological scales. Astronomical 
signals travel a long distance, L, and they thus correspond to an average 
over many oscillations. Instead of an oscillatory variation of photon number 
oc sin^(fcoscT/2), we would expect a constant reduction to a level given by 
the effective mixing angle in the dilute plasma between galaxies. 

The mixing of photons with hidden photons can occur in vacuum, but 
might be modihed in the presence of a plasma that gives the photons an ef¬ 
fective mass, whereas the WISP mass is essentially unchanged. This can give 
rise to matter oscillations reminiscent of the Mikheyev-Smirnov-WoIfenstein 
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(MSW) effect for neutrino oscillations |22UU221] . In particular, resonant con¬ 
version of photons into WISPs can occur within a plasma, i.e. both within 
astrophysical sources and during propagation of photons from the source to 
the observer. 

Once produced inside an astrophysical source, the weak coupling of WISPs 
to ordinary matter permits their escape without signihcant reabsorption, 
opening an additional channel of radiative cooling. The few neutrinos ob¬ 
served from the cooling phase of SN1987A have been used to place limits on 
the mass and coupling of axions by requiring that core-collapse supernovae 
not cool significantly faster than by neutrino emission alone |222j . 

Probing photon-WISP oscillations requires a detailed understanding of 
the photon-emission process. One of the best-understood radiation sources 
in the universe is the CMB. Apart from dipole distortion, its spectrum de¬ 
viates from a perfect blackbody by less than about 10“"^, and the amplitude 
of anisotropy is on the order of 10“®. The photons of the CMB have trav¬ 
elled a distance commensurate with the present-day Hubble radius. Any 
photon-WISP mixing at a level of would induce a spectral distortion 
or an anisotropy in conflict with CMB measurements. Since photon-ALP 
mixing requires the presence of a magnetic field, the absence of a significant 
modification of the CMB imposes an upper limit on the combination of the 
root mean square of the large-scale extragalactic magnetic-field strength and 
a parameter of the ALP model. 

In contrast to the case of ALPs, oscillations between photons and hidden 
photons do not require a magnetic field. The absence of CMB distortions 
thus directly constrains the vacuum mixing angle to be x 10“^ — 10“^ for 
hidden-photon masses of eV < my < 10“’^ eV |223] . 

The physics of non-thermal astrophysical sources is not well understood, 
resulting in a substantial uncertainty concerning the flux and spectrum of 
X-rays and y-rays. The observed photon spectra from these objects are 
often well represented by power laws, suggesting that the photons emissivities 
inside the sources are power laws as well. Then, photon-ALP mixing can 
induce steps in the spectra that may be detectable, and the absence of such 
spectral features can be used to place limits on ALP models. One of the free 
parameters of the theory is the Peccei-Quinn scale or axion decay constant 
fa |224] . Spectral modihcations on the keV to TeV scale are expected for ALP 
masses of m^ 10 ®eV and fa < 10^^ GeV |22511226] . somewhat depending 
on the magnetic-held strength. Exploiting these effects is complementary to, 
and potentially more sensitive than, other experimental avenues. 
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The currently best constraints have been derived from solar observations. 
Photons from the Sun would convert to ALPs in the solar magnetic held and 
then turn back into photons while passing through the held of an artihcial 
magnet in front of a telescope. The CERN Axion Solar Telescope (CAST) 
experiment provided the strongest constraint to date, /„ > 10^° GeV for 
ma < 0.02eV ^227\ . 

As photon-ALP mixing is energy-dependent, the resulting modihcation of 
the photon hux, or luminosity, will not be the same at diherent energies. One 
of the standard tools of astronomy is comparing the luminosity in diherent 
energy bands, which often reveals correlations that are specihc for a given 
class of objects. Photon-ALP mixing will impact the scatter of such corre¬ 
lations and induce non-Gaussianity in the distribution, an example of which 
was found in AGN data of optical/UV and X-ray emission |228] . Besides the 
obvious point that conventional astrophysics may lead to a non-Gaussian 
scatter distribution [229], substantial uncertainty in the interpretation arises 
from the question where exactly the conversion to ALPs occurs. If that hap¬ 
pens inside the AGN, one would deduce nia 10“^eV and fa — 3x 10® GeV. 
Most AGNs reside in galaxy clusters, and it is easy to see that photon-ALP 
conversion may also proceed in the /xG-level magnetic held in the intracluster 
medium. In that case the appropriate parameters would be ma 10“^^ eV 
and/, <10^0 GeV. 

Another possible signature for photon-ALP mixing is a modihcation of 
intergalactic 77 absorption [23011231] 1232] . Observations of AGNs at various 
redshifts suggest that the absorption of TeV-band 7 -rays by pair production 
with the extragalactic infrared background [233] is weaker than estimated. 
One possible explanation is the conversion of 7 -rays into ALPs at the source 
and reconversion into photons in the galactic magnetic held. Alternatively, 
conversion and reconversion may be caused by the much weaker extragalac¬ 
tic magnetic held, in which case one would deduce a completely diherent 
viable range for the parameters m, and /,. Besides, conventional astro¬ 
physics can oher explanations for weak 7 -ray absorption as well, for example 
if the sources actually produce cosmic rays of very high energy that along 
the entire line of sight produce photons through pion production in py col¬ 
lisions [2341 1235] 1236] . A requirement is a weak dehection of cosmic rays, 
implying a weak extragalactic magnetic held. This interpretation would be 
challenged if variability in the multi-TeV 7 -ray hux from distant AGN was 
observed. 
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5. LHC searches 

Due to their lack of electromagnetic interactions, non-baryonic nature and 
therefore (at most) weak Standard-Model interactions, dark-matter particles 
escape detection at colliders in a similar way as neutrinos. As a consequence, 
they produce a characteristic signal of missing energy. At hadron colliders 
such as the Large Hadron Collider (LHC) at CERN, only the missing trans¬ 
verse energy,can be observed, since the longitudinal momentum fractions 
of the incoming partons in the colliding hadrons are unknown. The missing 
Et is then determined from the recoiling observed objects such as jets, heavy 
quarks, photons and leptons. 

5.1. Model-dependent searches 

Assuming a particular dark-matter particle and its interactions with Standard- 
Model particles, it is possible to perform complete signal and background 
Monte Carlo simulations. Specihc kinematic cuts can then be dehned that 
largely eliminate the Standard-Model background and thus permit confronting 
the model expectations with experimental data. An excess in the kinematic 
distributions may then be attributed to the dark-matter particle in question. 

The background contribution can alternatively be determined directly from 
the data by extrapolating from a control region to the signal region. This 
procedure would, of course, not eliminate irreducible backgrounds from, e.g., 
neutrinos. 

If the dark-matter particle is the lightest and only stable member of a 
family of new states, it would be the only one having survived the evolu¬ 
tion of the universe. In that case, the typical LHC signal would consist of 
cascade decays of the heavier new states, thereby producing multi-parton or 
multi-lepton hnal states together with If the new spectrum includes 

strongly interacting particles, the corresponding hadronic cross sections are 
high, while the Standard-Model background rates decrease with the particle 
multiplicity in the hnal state. For this reason, searches for dark matter in 
cascade decays were the hrst to be performed at the LHC. 

An important example is the search for the lightest supersymmetric par¬ 
ticle (LSP), typically a neutralino in i?-parity conserving supersymmetry, in 
the pair production and subsequent cascade decay of squarks and gluinos 
to jets and two LSPs. Assuming a constrained Minimal Supersymmetric 
Standard Model (cMSSM), it is then possible to exclude regions in the pa¬ 
rameter space of mo and mi/ 2 , the common scalar and fermion cMSSM mass 
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Table 1: Operators coupling two dark-matter particles X to Standard-Model quarks (g) 
and gluons (g). The operator names beginning with D, C, and R apply to Dirac fermions, 
complex and real scalars, respectively. 


Name 

Operator 

Coefficient 




D1 

XXqq 

rriqlMf 




D2 

X-f^Xqq 

ifUq/Mf 

Name 

Operator 

Coefficient 

D3 

XXqj^q 

iniq/Mf 

Cl 

X^Xqq 

rriq/Mf 

D4 

X'j^Xq'j^q 

rriq/M^ 

C2 

XtXgySg 

imq/Mf 

D5 


l/Mf 

C3 

X^d^XqYq 

l/Mf 

D6 


l/Mf 

C4 

Xta^Xgy^ySg 

l/Mf 

D7 


l/Mf 

C5 

XtXG^.G'^" 


D8 


l/Mf 

C6 

XtXG^.G'^" 

ias/AMf 

D9 

Xa^^'^Xqa^yq 

l/Mf 

R1 

X^qq 

mq/2Mf 

DIO 

Xa^^'-f^Xqa^yq 

ijMl 

R2 

X^qj^q 

imq/2Mf 

Dll 


«,/4M3 

R3 

X^G^^G^" 

as/SMf 

D12 



R4 

X^G^^G^" 

ias/SMf 

D13 


ias/^Mf 




D14 


«s/4M3 





parameters that in general also determine the LSP mass |23711238] . 

An alternative is to search for direct production of dark-matter particles, 
which is observable only if it involves at least one visible particle in so- 
called mono-object events. This object can be either a jet, a photon, or an 
electroweak gauge or Higgs boson, and it is radiated either from the initial 
state or from an intermediate particle. The dark-matter signal then directly 
appears as an excess in the tail of its Et distribution. A typical example is 
the associated production of a jet with a pair of gravitinos, which are often 
the LSPs in gauge-mediated supersymmetry-breaking models |239j . 

5.S. Effective field theory approach 

A model-independent way to search for dark matter at the LHC is to 
analyze mono-object events within an effective held theory (EFT). Here, the 
coupling of dark matter to Standard-Model particles is parameterized with 
a set of non-renormalizable, higher-dimensional operators. The operators 
containing two dark-matter (fermionic or scalar) particles can be categorized 
and are shown in Tab. a BUI- As in Fermi’s four-fermion weak-interaction 
model, the effective couplings of mass dimension (—2) for these operators are 
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obtained by integrating out the propagator of the mediator with mass M and 
coupling strength g from j — M^) to g'^/M'^ =: 1 /M^ at low momentum 
transfer, p. Some of the coefficients in Tab. [T] are proportional to the quark 
mass niq and are thus more accessible in heavy-quark jets. 

In contrast to many explicit dark-matter models like supersymmetry, the 
EFT approach depends on only very few parameters, i.e. the scale of the 
interaction, M*, and the dark-matter mass, mx- Due to the crossing sym¬ 
metry of the underlying Feynman diagrams, it can also serve as a common 
standard for the comparison of LHC results with those from direct and in¬ 
direct searches cm- A necessary condition for the applicability of EFT is, 
however, that M* is indeed sufficiently large with respect to the momentum 
transfer |24211243] . The impact of this condition can be very different in di¬ 
rect detection and at the LHC. Second, if the mediator is light enough to be 
produced at the LHC, the collider analysis should be modihed. Conversely, 
if the dark-matter particle in question constitutes only a fraction of the local 
dark-matter density, the direct-detection limits change. Finally, the informa¬ 
tion obtained with EFT is of course less complete than that of a dedicated 
analysis. 


5.3. Experimental results from monoobject searches 

Recent monojet analyses by ATLAS |244j and CMS 


demonstrate 


that the LHC is quite competitive with direct dark-matter searches, in partic¬ 
ular for low masses. An overview of results is presented in Fig. focussing in 
particular on the operators Dl, D5, Dll, Cl and C5 for the spin-independent 
case and the models D8 and D9 for the spin-dependent case |244j . 

The monophoton analyses of ATLAS |246j and CMS |247j are somewhat 


less competitive than the monojet and also the monolepton results, as Fig. 16 


shows. The monolepton channel may furthermore involve different couplings 
of the dark-matter particle to up- and down-type quarks, and so the en¬ 
suing interference effects have been parameterized in the couplings with an 
additional factor ^ = ±1. 

In the future, measurements of correlations of visible particles, e.g. of 
the azimuthal angle between two jets, may permit more discriminative cuts, 
in particular between the signal and irreducible background coming from 
invisible Z decays to two neutrinos |248j . With search strategies adapted to 
the new pile-up and detector conditions, the LHC experiments should surpass 
the previous limits from center-of-mass energies of y/s = 7 and 8 TeV within 
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Figure 15: 90%-confidence upper limits obtained from a recent ATLAS monojet analysis at 
a center-of-mass energy of 8 TeV on the spin-independent (left) and spin-dependent (right) 
dark-matter/nucleon scattering cross section as a function of dark-matter mass mx for 
different operators. Results from direct-detection experiments for the spin-independent 
and spin-dependent cross section and the CMS results are sho'wn for comparison |244) . 


one year of data taking, reaching interaction scales M* of abont 2 TeV at 
-y/s =14 TeV with an integrated Inminosity of abont 100 fb“^ |249] . 

5.4- Simplified models 

An intermediate level of interpretation of LHC dark-matter searches are 
simplihed models. They provide microscopic descriptions of interactions be¬ 
tween dark matter and Standard-Model particles in a way that is still inde¬ 
pendent of (nltraviolet) complete extensions of the Standard Model. They 
can be classihed into fonr basic types of messenger helds: scalar, psendo- 
scalar, vector, or axial-vector. In this way, they are still characterizable by 
a relatively small number of four to hve parameters, including the mediator 
mass and width, the dark-matter mass and one or two effective couplings 
[250]. 

The ATLAS |251j and CMS [247] collaborations have searched at ^/s = 8 
TeV for the production of single Z bosons and single photons in association 
with missing Ex, respectively. Z bosons are identihed by their leptonic (i.e. 
oppositely charged electron or muon) decay products. Since no excess above 
the Standard Model predictions was observed, limits were hrst set on the mass 
scale of the contact interaction, M* (see above), but also on the coupling and 
mediator mass in simplihed models. 
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Figure 16: 90%-confidence upper limits obtained from a recent CMS monophoton analysis 
at center-of-mass energies of 7 and 8 TeV on the X-nucleon cross section as a function of the 
dark-matter particle mass, Mx, for spin-independent (left) and spin-dependent couplings 
(right) [247j . Also shown are limits from CMS using monojet and monolepton signatures, 
including results for maximized interference parameter, ^ = ±1, as well as several published 
limits from direct-detection experiments. The solid/hatched yellow regions show the 68% 
and 95%-confidence contours respectively for a possible signal in CDMS. 


In the ATLAS analysis, only scalar mediators were considered. Figure [T7| 
(left) shows as a function of coupling, /, the permitted regions in the param¬ 
eter plane spanned by the dark-matter mass, mx, and the scalar-mediator 
mass, m^. Also included are limits in the mx-m^j plane obtained by impos¬ 
ing the observed dark matter relic abundance (thick black line), which in the 
analysis was possible for dark matter and mediator masses up to 325 GeV 
and 1200 GeV, respectively. 

In the GMS analysis, only vector mediators were considered. Figure 17 


(right) shows, as a function of the mediator mass, limits on the mediator 
mass, M, divided by the geometric average of the coupling to dark matter, 
gx, and that to quarks, gq. 

For Run II of the LHG, simplihed models will play a prominent role. 
Experimental signatures of the types and characteristics of mediators, e.g., 
their couplings to dark matter and Standard-Model particles, will be used as 
building blocks for a common model implementation, benchmark definitions, 
and assessments of the validity of the EFT approach. For the monojet sig¬ 
nature, for example, resonant s-channel vector and non-resonant f-channel 
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Figure 17: Left: 95%-confidence upper limits obtained from a recent ATLAS mono-^- 
boson analysis at a center-of-mass energy of 8 TeV on the coupling constant, /, of the 
scalar-mediator theory as a function of mx and the mediator mass, m,,. The cross- 
hatching indicates the theoretically accessible region outside the range covered by this 
analysis. The white region is phase space beyond the model’s validity. In the excluded 
region in the upper left-hand corner, demarcated by the black line, the lower limit on 
/ from the relic-abundance calculations is larger than the upper limit measured in this 
analysis |251j . Right: Observed limits obtained from a recent CMS monophoton analysis 
at a center-of-mass energy of 8 TeV on the dark-matter mediator mass divided by coupling, 
M/y/gxgq, as a function of the mediator mass, M, assuming vector interactions, for dark- 
matter particle masses of 50 GeV and 500 GeV. The width, T, of the mediator is varied 
between M/Sir and M/3. The dashed lines show contours of constant coupling |247j . 


scalar mediators coupling to quarks will be considered together with scalar 
mediators coupling to gluons through heavy (in particular top) quark loops 

[219]. 

A rationale for this choice are models such as the inert-Higgs doublet 
|252] . which is motivated by the recent Higgs discovery and the possibility 
that dark matter primarily couples to Higgs bosons. The advantage of mini¬ 
mal models lies in their small parameter space and complete testability, not 
only at the LHC, but also in direct-detection experiments [253]. Signatures 
including a top quark will also be considered on account of its large coupling 
to Higgs bosons and its important role in electroweak symmetry breaking. 

6. Towards global analyses 

In this article, we have reviewed the observational evidence for the pres¬ 
ence of dark matter in our universe, covering a vast range of astrophysical 
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scales, and possible candidates from elementary particle physics together 
with their theoretical relations to various proposed extensions of the Stan¬ 
dard Model. We have then summarized the three main roads followed in 
dark-matter searches: direct detection, based on interactions of the local 
dark matter in dedicated laboratory experiments, indirect detection, based 
on the search for decay and annihilation products of dark matter on astro- 
physical and cosmological scales with multiple messengers, and searches at 
the LHC. Special emphasis was placed on the complementarity of these three 
detection modes. All three have specific advantages in that they can provide 
critical information on specihc particles that might constitute dark matter. 
Finding and measuring the properties of GeV- to TeV-scale new particles at 
the LHC is a unique channel for exploring the interactions with Standard- 
Model particles. Scattering of dark-matter particles off various materials in 
direct-detection experiments can provide insight into the local density and 
velocity distribution of dark matter. Indirect detection through astrophysical 
measurements permits tracing the distribution of dark matter in the sky. 

Finding new particles that constitute dark matter would be a major 
breakthrough in physics. As extraordinary new hndings require extraordi¬ 
nary evidence, the hurdles are high. Various unusual experimental results of 
the last decade have been interpreted in terms of dark matter, but all of them 
could also be the result of a misunderstood background or other effects. The 
problem is probably most obvious for the indirect detection in astrophysical 
messengers, where conventional astrophysical processes can lead to surpris¬ 
ing observational features. It is therefore mandatory that all three strategies 
of dark-matter detection be continued. There is complementarity in the de¬ 
tection modes, but there is also overlap in the sense that they sample similar 
parts of the parameter space with independent methods. The best option 
for the future appears to be a global approach, in which the three detection 
modes are jointly pursued and the parameter space of effective operators, 
classes of simplihed models, and eventually specific models are tested. Com¬ 
munication between the now often separately operating research initiatives 
will be a key for success. Finding and identifying dark matter requires thus 
a world-wide concerted effort. 

We have also briefly discussed alternatives to dark matter, that may ex¬ 
plain part of the observational indications of excess gravitational action. Al¬ 
though the evidence for dark matter composed of individual particles is huge, 
in science one is well advised to keep an open eye to alternative or additional 
physical processes that may be at play. 
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